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ABSTRACT 


\ 

\ 

\  / 

;i 

Basic  theories  for  polarization  utilization  in  radar  target 
reconstruction  are  presented  and  a  general  literature  review  with 
many  pertinent  references  is  given. 

The  mathematical  descriptions  of  polarization  in  terms  of  the 
Poincare  polarization  sphere  are  introduced  and  the  relationships 
existing  among  the  radar  scattering  matrix  [S],  the  Stokes 
reflection  matrix  [M] ,  the  modified  Mueller  matrix  [Mm],  and  the 
coordinates  of  the  related  co-polarization  and  cross-polarization 
nulls  on  the  Poincare  sphere  are  derived. 

It  is  shown  that  a  scattering  phenomenon  can  be  uniquely 
expressed  given  the  elements  of  either  one  of  [S],  [M] ,  [Mm]  or 

the  coordinates  of  the  optimal  polarizations,  i.e.  unique 
inversion  relations  among  the  four  equivalent  representations 
exist  which  is  relevant  to  target  polarization  synthesis. 


The  developed  theories  are  verified  by  computer  computation 
using  measurement  data  and/or  model  scattering  data  as  inputs. 
The  compucer  programs  are  listed  and  examples  of  our  optimal 
polarization  analysis  are  presented  for  the  monostatic,  relative 
phase  case.  Single  perfectly  conducting  target  shapes  and  some 
sea  clutter  testing  models  witn  and  without  target  were  chosen; 
and  our  studies  demonstrate  clearly  that  the  optimal  polarization 
concept  introduced  by  KENNAUGH  is  very  useful  in  radar  target 
analysis  as  will  be  further  pursued  in  ether  forthcoming  reports. 
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CHAPTER  ONE 


INTRODUCTION 

In  recent  yearst  there  has  been  a  rapidly  expanding  volume  of 
research  from  both  a  theoretical  and  experimental  point  of  view, 
directed  towards  the  determination  of  the  characteristic 
properties  of  radar  targets  through  the  use  of  polarization.  The 
fact  that  makes  this  type  of  investigation  possible  is  that  the 
scattering  properties  of  radar  targets  are  dependent  on  the 
polarization  of  the  incident  radiation.  This  dependence  which 
manifests  itself  as  a  depolarization  of  the  scattered  wave,  is  a 
function  of  the  structure  and  geometry  of  the  scatterer  thus 
beirg  characteristic  for  a  particular  target.  It  has  been 
demonstrated  that  a  radar  target  acts  as  a  polarization 
transformer.  This  transformation  was  expressed  (Sinclair  1948)  as 
a  matrix  [S]  which  could  be  Incorporated  into  the  radar  range 
equation.  Kennaugh  (1950)  gave  a  geometrical  meaning  to  the 
transformation  by  representing  the  power  received  by  a  radar  on 
the  Poincare  sphere.  In  his  series  of  reports  (1949-1954), 
Kennaugh  also  demonstrated  that  there  exist  two  radar 
polarization  states  for  which  the  radar  receives  zero  signal  from 
the  target!  which  are  known  as  null  polarizations.  As  it  will  be 
shown  later  the  null  polarizations  (co-and  cross-polarization 
nulls)  along  with  the  polarization  orientation  invariant,  i.s. 
the  Span  of  the  scattering  matrix  [S]  (span{ [S] }-£^  I ^ 

| Sij | 1  .where  the  are  the  elements  of  the  matrix  [S]), 

can  be  used  to  describe  in  their  entirety  the  characteristic 
properties  of  a  target  at  one  frequency  and  for  one  aspect. 
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It  has  been  established  experimentally  (Daley  1978-79, 
Weisbrod  and  Morgan  1979)  that  the  null  polarizations  can  be  used 
in  order  to  discriminate  targets  against  scattering  clutter, 
which  is  of  great  importance.  This  meets  greatly  the  existing 
tremendous  need  for  improved  clutter  rejection  methods  in  order 
to  detect  accurately  small  surfaces  which  is  also  of  great 
significance  to  military,  geophysical  and  environmental  remote 
sensing.  It  has  been  shown  (Weisbrod  et  al,  1979)  in  the  case  of 
sea  clutter  that  its  non-random  behaviour  manifests  itself  as  a 
characteristic  clustering  of  co-and  cross-polarized  nulls  as 
plotted  on  the  Poincare  sphere.  This  clustering  was  noticeably 
disturbed  with  the  presence  of  a  target.  This  phenomenon  could 
lead  to  zero  false  alarm  rate  discriminants  with  the  use  of 
theoretical  models  extending  existing  clutter  statistics, 
sensitive  to  the  changes  in  the  clustering  of  co-pol  and  cross- 
pol  nulls. 

In  view  of  the  fact  that  the  co-and  cross -polarization  nulls 
along  with  the  polarization  transformation  invariant  (span  of 
matrix  [S]  )  are  characteristic  of  a  particular  typ>  of  target 
for  one  aspect  and  at  one  frequency,  these  quantities  should  be 
given  directly  from  the  measurement  data  and  recorded  on  the 
polarization  sphere  and/or  its  associated  polarization  maps. 

So  far,  the  polarization  utilization  in  radar  target  related 
phenomena  has  been  accomplished  through  the  use  of  the  radar 
scattering  matrix  [S]  and  the  optimal  polarization  pairs  of  the 
co-and  cross-polarization  nulls  have  been  expressed  in  term  of 
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the  [S]  matrix  components.  The  use  of  the  Stokes  reflection 
matrix  or  Mueller  matrix  (Mueller,  1948)  has  been  overlooked  even 
though,  it  contains  as  much  information  as  the  scattering  matrix 
and  furthermore  it  is  easier  to  obtain,  since  it  involves  only 
amplitude  (power)  measurements  and  not  phase  determination  as  in 
the  case  for  the  matrix  [S],  There  is  also  convincing  evidence 
(Leader  1978,  Boerner  1979)  that  the  Mueller  matrix  elements 
behave  in  a  manner  characteristic  of  the  material  properties  of 
the  specific  type  of  clutter  they  represent.  Though  the 
interpretation  of  this  behaviour  is  in  need  of  further 
investigation,  one  can  safely  assume  that  the  information 
provided  b  ;he  matrix  [M]  can  be  very  useful  in  establishing 
target  and  particularly  clutter  characterisitica  especially  for 
the  incoherent  scattering  case.  Due  to  the  "additivity"  property 
(Chandrasekhar  1950,  Ishimaru  1978)  of  the  Stokes  parameters 
(Stokes,  1852)  of  independent  waves,  the  independent  incoherency 
properties  of  clutter  are  explicitly  contained  in  the  matrix  [M] 
while  they  are  only  implicitly  inherent  in  the  [S]  matrix.  In 
addition,  most  of  the  experimental  statistical  clutter 
distributions  are  given  ,  in  terms  of  the  components  of  [M] . 
Therefore,  in  order  to  interpret  scatter  characteristic 
polarization  properties  in  terms  of  the  associated  co-and  cross¬ 
polarization  null  statistics  on  the  Poincare  sphere,  it  would  be 
highly  practical  to  express  these  optimal  polarizationa  in  terms 
of  the  Mueller  matrix. 

In  view  of  the  above,  we  have  obtained  here  the  scattering 
matrix  elements  from  the  Mueller  matrix  (both  the  matrix  [M]  and 
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its  modified  form  [Mmj  for  the  bistatic  and  monostatic  cases). 
Thus  the  problem  of  obtaining  the  scattering  matrix  with  relative 
phases  from  the  knowledge  of  the  Mueller  matrix  is  feasible.  The 
expressions  related  to  the  problem  were  subsequently  tested 
against  real  measurements  for  the  case  of  sea  clutter  provided  to 
us  by  Mr.  J.  Daley  of  the  Naval  Research  Laboratory.  In  addition, 
we  have  shown  that  the  elements  of  [S]  ,[M]  and/or  [Mm]  can  be 
generated  from  the  knowledge  of  the  radius  of  the  Poincaro  sphere 
and  from  the  spherical  coordinates  of  either  both  co-polarization 
(COPOL)  or  one  COPOL  plus  one  cross -polarization(XPOL)  null.  Thjs 
inverse  relationship  is  important  as  it  will  assist  us  greatly  in 
the  target  polarization  synthesis  problem. 

The  text  of  the  present  report  is  comprised  of  five  Chapters. 
In  Chapter  Two  a  general  review  of  the  existing  theoretical  and 
experimental  work  in  radar  target  polarization  related  phenomena 
is  given.  Chapter  Three,  covers  the  theoretical  principles  on 
which  the  Basic  Polarization  Descriptors,  the  Scattering  matrix 
[S],  its  Transformation  Invariants  and  Mueller  Matrix  [M]  are 
based.  In  Section  3.6.2  of  the  same  Chapter  it  is  shown  how  the 
amplitudes  and  relative  phases  of  the  elements  of  the  scattering 
matrix  [S]  are  derived,  given  the  Mueller  matrix  [M]  or  its 
modified  form  [Mm] .  Also,  in  the  same  Chapter  the  relationship 
between  [M]  and  [Mm]  is  given.  Similarly,  we  show  how  the 
elements  of  [S]  ,  [M]  and/or  [Mm]  can  be  regenerated  from  the 
knowledge  of  the  spherical  coordinates  of  the  two  COPOL  nulls  or 
one  COPOL  plus  one  XPOL  null.  In  Chapter  Four  the  expressions 
obtained  in  Chapter  Three,  were  tested  satisfactorily  against 
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the  data  available  to  us.  Also  for  the  sake  of  completeness,  we 
have  calculated  in  Chapter  Four  the  Mueller  matrix  elements  (for 
both  [M]  and  [Mm]  matrices  in  the  monostatic  case)  from  the 
scattering  matrix  elements  generated  for  various  target  shapes 
and  sea  clutter. In  addition,  the  optimal  polarization  pairs  of 
the  co-and  cross -polarized  nulls  were  calculated  and  their 
coordinates  on  the  Poincare  sphere  were  found.  Finally  the 
summary  of  our  results,  conclusions  and  recommendations  are  given 
in  Chapter  Five.  Proofs  of  some  pertinent  identities  and 
computer  programs  for  calculating  relevant  parameters  are  given 
in  the  appendices . 
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CHAPTER  TWO 
GENERAL  REVIEW 

The  utilization  of  polarization  in  radar  target  and  clutter 
studies  has  been  the  subject  of  rather  extensive  theoretical  and 
experimental  research  efforts  in  recent  years . 

First  Sinclair (1948)  showed  how  a  radar  target  can  act  as  a 
polarization  transformer  and  he  was  able  to  express  this  property 
with  the  use  of  a  matrix  which  was  incorporated  in  the  radar 
range  equation.  There  followed  a  series  of  papers  by  Booker, 
Rumsey,  Deschamps,  Kales  and  Bohnert(195 1)  on  polarization,  as  it 
is  related  to  radar  antennas,  which  constituted  the  basis  for 
future  research  on  the  subject.  Huynen  et  al(1953)  initiated  from 
then  on  studies  on  radar  returns  by  investigating  the  effects  of 
polarization  on  radar  scattering  by  ground  targets  and 
precipitation.  Among  the  pioneers  in  the  field  of  theoretical 
work,  on  radar  target  scattering  were  Kennaugh( 1949- 1954)  and 
Gent(1954).  Kennaugh  was  also  the  first  worker  to  introduce 
concepts  of  such  potential  impact  as  the  optimal  polarizations  of 
a  radar  target(1949) .  Graves  (1956)  introduced  the  polarization 
power  matrix  which  gives  the  total  power  back-scattered  from  the 
target  for  any  transmitted  polarization.  Based  on  Kennaugh' s 
work,  Copeland(1960)  gave  a  classification  of  the  single  radar 
target  by  measuring  the  received  complex  voltages  using  rotating 
linearly  polarized  antennas.  Studies  on  polarization 
characteristics  in  the  scattering  from  symmetric  radar  targets 
were  reported  by  Crispin( 1961) ,  Bechtel  .and  Ross (1962).  and 
Huynen(1962,  1963).  A  summary  of  updating  radar  measurements  was 


presented  by  Huynen,  Landry,  Webb  and  Allen  at  the  Radar 
Reflectivity  Measurements  Symposium  in  1964.  In  a  Special  IEEE 
Proceedings  Issue  on  radar  reflectivity  (August  1965),  the  theory 
and  measurement  techniques  for  target  scattering  matrices 
(asymmetric  objects  included)  were  discussed  by 
Lowenschuss (1965) ,  Huynen(1965) ,  Copeland(1960) ,  Kuhl  and 
Covelli(1965) . 

A  vast  amount  of  research  has  also  been  performed  on  time- 
varying  distributed  targets  which  is  independent  of  the 
previously  cited  work  on  single  radar  targets.  In  his  analysis  of 
single  targets,  Gent(1954)  considered  also  ensembles  of 
distributions  of  single  targets.  Ament (1960)  studied  the  problem 
of  whether  reciprocity  holds  for  rough  surface  scattering  and 
Ko(1962)  gave  an  introduction  with  application  to  partially 
polarized  scattering.  Several  statistical  models  for  terrain  are 
given  by  Spetner  and  Katz(196G).  Studies  of  scattering  from  rough 
surfaces  using  scalar  theory,  which  will  remain  classical,  were 
treated  by  Beckmann  and  Spizzichino(1963)  and  updated  by 
Beckmann(1965) .  Various  rough  surface  models  related  to  both 
theory  and  measurement  have  been  studied  by  Beckmann (19 65 ) , 
Parks(1964)  and  Renau  and  Collinson(l965) .  Fung(March  1966,  July 
1966,1967)  treated  rough  surface  scattering  using  vector  thoery 
which  also  considered  depolarization  of  electromagnetic  waves. 
Using  high-frequency  asymptotic  theories  based  on  the  Kirchhoff 
approximation,  Hagrors(1967)  studied  the  depolarization  of  Lunar 
Radar  Echoes  and  Stogryn( 1967)  worked  on  electromagnetic 
scattering  from  rough, finitely  conducting  surfaces.  Van  de 


PAGE  8 


Hulst(1957)  gave  an  introduction  to  high-frequency  scattering. 
Worth  mentioning  here,  are  also  the  overall  reviews  of 
contributions  to  wave  scattering  from  rough  surfaces  and  various; 
kinds  of  clutter,  given  in  the  Transactions  of  the  IEEE  Special 
Issue  on  Partial  Coherence(1967) ,  in  "Radar  Astronomy"  edited  by 
Evans  and  Hagfors(1968) ,  in  Bowman  et  al(1969),  Crispin  and 
Siegel(1968) ,  in  Ruck  et  al(1970)  and  in  Beckmann's  book  on  the 
depolarization  of  electromagnetic  waves(1968).  Huynen(1970) 
developed  a  phenomenological  theory,  applicable  to  all  radar 
targets,  according  to  which  the  radar  target  appears  as  an  object 
for  investigation  through  the  process  of  the  radar  illumination. 
Thus  one  can  avoid  going  into  constructions  of  specific 
statistical  and  geometrical  models,  as  it  has  been  the  case  in 
most  of  the  literature  on  the  subject.  The  mathematical  framework 
for  the  theory  is  given  in  terms  of  the  target  polarization 
scattering  matrix,  which  has  been  used  extensively  in  the  Russian 
literature(Kanareykin  et  al,  1966,  1968;  Stead,  1967; 
Zhivotovskiy ,  1973;  Potekhin  et  al,  1969;  Kozolov,  1979,  and 
Basalov  et  al,  1973). 

Due  to  its  importance,  several  efforts  have  been  made  recently 
to  use  the  polarization  sensitivity  of  radar  targets  towards  a 
classification  of  radar  targets.  Thus  the  classification  of  radar 
targets  suggested  by  Copeland  has  been  recently  investigated 
experimentally  by  Steinbach(1973 ,  1976).  Studies  on  some  methods 
for  radar  target  identification  are  given  by  Von  Schlachta(1977) , 
Crom(1973),  and  Jeske(1976).  There  has  also  been  or.  extensive 
amount  of  work  to  be  found  in  the  recent  literature  on 
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precipitation  scatter  (Schneider  and  Williams,  1976;  McCormick 
and  Hendry  1976;  Barge  and  Humphries,  1979,  1980;  Hall  et  al, 
1980;  Bringi  et  al,  1979),  on  scattering  from  irregular 
surfaces (Bass  and  Fuks ,  1975)  and  sea  clutter (Valenzuela  1968, 
1978).  A  detailed  review  on  radar  reflectivity  of  land  and  sea 
is  given  in  the  textbook  of  Long(1975)  and  also  was  discussed  in 
Stiles  and  Holtzmann,  1979. 

Of  particular  interest  to  the  present  work  is  the  recently 
developed  literature  on  studies  of  the  co-and  cross -polarization 
nulls  of  radar  targets  which  promise  to  lead  to  very  effective 
radar  target  versus  clutter  discrimination  techniques.  Recently, 
studies  made  in  the  Naval  Research  Laboratory (Daley  1978,  1979; 
Weisbrod  and  Morgan  1979)  showed  that,  in  the  case  of  sea  clutter 
the  co-and  cross-polarization  nulls  of  clutter  exhibit  a  non- 
random  stable  distribution  when  mapped  on  the  Poincare  sphere. 
This  distribution  was  disturbed  in  the  presence  of  a  target,  a 
phenomenon  that  may  lead  to  effective  target  discrimination  when 
fully  investigated.  The  importance  of  the  co-and  cross¬ 
polarization  nulls  in  radar  target  versus  clutter  discrimination 
has  been  emphasized  particularly  by  Poelman(1971 ,  1974,  1976, 
1977).  Recently  Ioannidis  and  Hammers ( 1979) ,  Rosien  et  al(1979) 
have  suggested  schemes  for  radar  identification  in  clutter,  based 
on  the  characteristic  properties  of  the  co-polarization  nulls. 
We  also  wish  to  refer  to  the  forthcoming  IEEE  Trans.  AP  29(2), 
March  1981,  Special  Issue  on  "inverse  Methods  in 
Electromagnetics",  in  which  several  papers  on  polarization 
correction  are  presented. 
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Target  versus  clutter  discrimination  for  targets  embedded  in 
terrain  clutter  has  been  improved  by  utilizing  also 
depolarization  effects (Gent  et  al,  1963;  Weiss,  1967;  Egan  et  al, 
1967;  Dainty,  1975).  It  was  shown  by  Gough  and  Boemer(1979) , 
Egan  et  al(1967)  that  the  correlation  functions  of  linearly  and 
circularly  polarized  components  allow  discrimination  between 
metal  and  dielectric  random  scatterers.  In  their  studies  Gough 
and  Boerner (1979)  are  using  the  additive  properties  of  the 
Mueller  matrix  in  order  to  analyse  the  interaction  of  coherent 
target  signal  with  various  incoherent  clutter  components.  The 
usefulness  of  the  Mueller  matrix(Mueller  1948,  McCormick  1950) 
has  also  been  recognized  by  Hagfors (1967) ,  Leader (1978)  and 
others . 

Concluding  we  cite  here  other  textbooks  and  papers  related  to 
radar  target  scattering  as  well  as  to  the  theory  on  partial 
coherence.  Thus  we  distinguish  here  the  works  by  Kraus(1966), 
Nathanson(1969) ,  Thiel(1970),  Meyer  and  Mayer(1973),  Boerner  (The 
State  of  the  Art  Review,  1978),  Born  and  Wolf (1966), 
Ishimaru(1978) ,  Strohbehn(1978) ;  in  the  Russian  literature  the 
works  by  Kobzarev(1969) ,  Shirman(1970) ,  Gorshkov(1974) , 
Kozlov(1979)  and  Zhivotovskiy(1976) .  We  note  here  that  in 
optics  major  contributions  were  made  by  R.  Clark 
Jones ( 1941-1956)  as  was  documented  in  Shurclif f ( 1962) ,  Hecht  and 
Zajac(1976)  and  Roots(1978). 
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CHAPTER  THREE 
THEORY 

3.1  Introduction  : 

Chapter  Three  discusses  iji  general  the  theoretical  principles 
related  to  the  utilization  of  polarization  in  radar  target 
scattering  phenomena.  First,  in  Section  3.2,  the  basic 

descriptors  of  the  polarization  of  monochromatic  plane  waves  are 
given.  In  Sections  3.3  and  3.4,  the  theory  of  the  radar  cross- 
section  scattering  matrix  is  discussed  and  its  transformation 
invariants  as  well  as  its  representation  on  the  Poincare  sphere 
with  the  use  of  the  co-and  cross-polarization  nulls  are 
considered.  In  Section  3.5,  the  Stokes  or  Mueller  matrix  is 
reviewed  and  a  method  of  reconstructing  the  scattering  matrix, 
given  its  associated  Mueller  matrix  is  introduced  in  Section  3.6. 
In  the  same  Section,  the  relationship  between  the  Mueller  and 
modified  Mueller  matrices  is  given.  Finally  in  Section  3.7,  the 
elements  of  the  scattering  and  the  Mueller  matrices  are 
regenerated  from  the  knowledge  of  the  spherical  coordinates  of 
either  of  the  two  co-polarlzation(COPOL)  nulls  or  one  COPOL  plus 
one  cross-polarization(XPOL)  null. 

3.2  Basic  Polarization  Descriptors  : 

3.2.1  The  Polarization  Vector 

The  polarization  of  an  electromagnetic  wave  describes  the 
orientation  of  the  field  vectors  at  a  given  point  during  one 
period  of  oscillation.  In  the  present  treatment,  polarization  is 
referred  to  the  orientation  of  the  electric  field  vector  E  only, 


since  we  are  concerned  with  the  far-field  scattered  radiation  and 
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the  magnetic  field  vector  H  is  perpendicular  to  E  in  direction 
and  proportional  to  that  of  E  in  magnitude.  The  direction  along 
which  a  wave  propagates  is  given  by  the  propagation  vector  k.  In 
an  isotropic  medium,  the  plane  containing  E  and  H  and  which  is 
perpendicular  to  the  direction  of  propagation,  is  called 
polarization  plane  [IEEE  Standard  149,  1979].  In  the  case  where 
the  polarization  of  a  wave  is  the  same  at  every  point  in  space, 
the  wave  is  said  to  be  linearly  or  plane  polarized.  We  consider 
here  monochromatic  plane  waves,  i.e.,  plane  waves  at  a  single 
angular  frequency  u,  propagating  along  the  z-axis  of  a 
rectangular  coordinate  system  xyz,  where  the  xy  plane  is  a 
reference  plane  such  as  the  mean  surface  of  the  earth.  The 
electric  field  vector  of  a  time-harmonic  plane  wave  at  a 
position  r  (x,  y,  z)  and  time  t  is  given  by 


ECi.t)  =  =  he^  Cut-kz) 

where  k  is  the  propagation  vector  with  magnitude  k  =  2tt/X,  X  is 
the  wavelength  of  the  wave  in  free  space,  h  is  the  complex 
electric  field  phasor,  known  as  the  complex  polarization  vector. 
In  radar  propagation  the  vector  h  =  h”^,  with  complex  magnitude  h 
and  direction  specified  by  the  unit  vector  ia,  may  be  decomposed 
along  the  two  orthogonal  directions  x  and  y  represented  by  the 
unit  vectors  h^  and  h^,  in  the  following  manner  : 


(3.1) 


h  =  h  ^  =  h  /hs  fh'h' 
-  x  x  y  y 


(3.2) 


! 


( 
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where"  a  ,  a  are  the  polarization  vector  component 
x  y 

magnitudes,  6  ,  5  are  their  phases  and  6=6  -  6 

x  y  y  x 

is  their  phase  difference.  In  the  antenna  language,  the 

magnitude  of  the  polarization  vector  |  h  | 2  =  a2  =  a  2  + 

X 


ay2  is  a  measure  of  antenna  radiation  efficiency  and  6  is 


the  phase  difference  between  the  x  and  y  channels  of  the 
antenna  (Kraus,  1966). 


3.2.2  The  Polarization  Ellipse 

According  to  (3.1)  and  (3.3),  the  electric  field  vector  E 
consists  of  two  components  along  the  x  and  y  axes  which  are 
given  by  : 


E  -  R.{h  .Kut-kz), 

X  X 

=  a  cos(e  +  6  ) 

X  X 

E  =  Re{h 

y  y 


=  Re{a  e^x  e^  ^ } 
x 

=  Re{ayej6y  ej(Wt_kz)} 


(3.4a) 


(3.4b) 


or 

A  ^  A  A 

E=Eh+Eh  =  a  cos(e+6  )h  +a  cos(e+6  )h  (3.4c) 

•—  XX  y  y  X  X  y  y '  y  '  ' 

where  e  =  wt  -  kz,  Re{»)  stands  for  the  real  part  of  a 
complex  number  {•). 


The  curve  that  the  electric  field  vector  E  describes  at  a 
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typical  point  in  space,  is  the  locus  of  the  points  whose 
coordinates  are  given  by  (3.4).  The  nature  of  this  curve  is 
found  after  eliminating  e  between  Eqs.(3.4)  and  after  some 
algebraic  manipulation,  to  be  given  by  : 


which  is  the  equation  of  an  ellipse  known  as  the 

polarization  ellipse.  The  ellipse  is  inscribed  into  a 

rectangle  whose  sides  are  parallel  to  the  coordinate  axes 

and  whose  lengths  are  2a  and  2a  (Fig.  3.1). 

x  y 


In  Eq.  (3.4)  if  a^  =  0  ,  the  wave  is  linearly  polarized 


in  the  y-direction.  If  ay  =  0*  the  wave  is  linearly 

polarized  in  the  x-direction.  Following  Kraus (1966) ,  if  5  = 

% 

6  -6=0  and  a  =  a  ,  the  wave  is  also  linearly 

y  x  x  y 

polarized  but  in  a  plane  at  angle  of  45°  with  respect  to  the 
x-axis.  A  further  special  case  of  interest  occurs  when  a^ 


=  a^.  and  5  =  ±90°.  The  resulting  wave  is  circularly 

polarized.  When  5  =  +90°,  the  wave  is  said  to  be  left 

circularly  polarized  and,  when  5  =  -90°,  it  is  said  to  be 


right  circularly  polarized.  Thus,  from  (3.4)  we  have  for  6  = 

+90°,  at  5  =  0,  z  =  0  and  t  =  0,  that  E  =0  and  E  = 
y  x  y 


a  as  in  Fig.  3.2a.  Under  the  same  conditions  but  at  a 

y 

later  time  such  that  wt  =  90°,  E  =  0  and  E  =  a  ,  as 

y  xx 


shown  in  Fig.  3.2b.  The  rotation  of  the  electric-field 


PAGE  16 


vector  Is  thus  clockwise  with  the  wave  approaching. 
According  to  the  IEEE  standards (19 79)  this  sense  of  rotation 
is  definded  as  left  circular  polarization.  According  to  the 
older  usage  of  classical  physics  (Born  and  Wolf),  this  sense 
of  rotation  (clockwise  with  wave  approaching)  is  definued  as 
right  circular  polarization,  or  opposite  to  the  IEEE 
definition  (Kraus  1966). 

If  the  wave  is  viewed  receding  (from  negative  z  axis  in 
Fig. 3-1),  the  electric  vector  appears  to  rotate  in  the 
opposite  direction.  Hence,  clockwise  rotation  with  the  wave 
receding  is  the  same  as  counterclockwise  rotation  with  the 
wave  approaching.  In  the  following  the  IEEE  definition  will 
be  used,  since  it  could  also  be  defined  (without  reference 
to  the  wave  direction)  by  means  of  helical-beam  antennas 
(Kraus,  1950).  Thus,  a  right  handed  helical-beam  antenna 
radiates  or  receives  right  circular  (IEEE)  polarization.  A 
right-handed  screw,  is  right  handed  regardless  of  the 
position  from  which  the  helix  is  viewed.  There  is  no 
possibility  here  of  ambiguity.  In  general  Eq.  (3.5) 
represents  a  left-handed  elliptical  polarized  (ep)  wave  if 
sin  6  >  0  and  right-handed  if  sin  6  <  0. 

We  now  seek  ultimately  to  define  the  polarization  vector 
in  terms  of  the  geometric  descriptors  of  the  ellipse  of  the 
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elliptical  polarized  (ep)  plane  wave  it  represents.  The 
reason  for  this,  is  that,  Eq.  (3.3)  is  inconvenient  when  it 
comes  to  dealing  with  the  proporties  of  radar  targets  since 
h  is  being  associated  with  &  fixed  (xyz)  antenna  coordinate 
system  of  which  the  radar  targets  are  independent. 

Hie  geometric  descriptors  of  an  ellipse  are  : 

1)  its  size  given  by  the  magnitude  a  =  |  h  |  of  the 

polarization  vector, 

2)  the  orientation  of  the  ellipse  with  respect  to  the  x-axis 
given  by  the  orientation  angle  0(0  £  0  £  tt ) , 

3)  the  ellipticity  angle  t,  (-it/4  S  t  £  ir/4)  such  that  cot  t 
is  given  by  the  ratio  of  the  semimajor  to  semirainor  axes  of 
the  ellipse  and  it  represents  the  axial  ratio  (AR) ,  where 
ISA 8S«,  and 

4)  the  sense  in  which  the  ellipse  is  being  traversed.  It 
will  be  shown  later  that  the  sense  can  be  given  by  the  sign 
of  t.  Though,  it  seems  natural  to  define  the  sense  as  right- 
handed  or  left-handed  according  to  whether  the  rotation  of 
the  electric  vector  E  and  the  direction  of  propagation  from 
a  right-handed  or  left-handed  screw  as  explained  before. 

We  now  choose  a  coordinate  system  (x',y')  such  that  the 
ellipse  has  an  orientation  angle  0  =  0  in  these  coordinates. 
According  to  Fig.  3.1,  the  polarization  vector  h  in  this 
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case  can  be  given  by  : 


h(a,T,o) 


acost  jo 
e 

j  as  int 


(3.6) 


where  a  is  called  the  ''absolute  phase"  of  the  antenna,  and 
defines  the  phase  reference  of  the  antenna  at  time  t  =  0 . 
From  (3.6),  we  notice  that,  when  the  sign  of  t  changes,  the 
direction  of  sense  of  polarization  changes,  being  left- 
sensed  for  positive  values  of  t  and  right-sensed  for  the 
negative  ones. 


The  general  expression  of  the  polarization  vector  for  an 
elliptically  polarized  wave  with  an  ellipse  of  orientation  <j> 
is  obtained  from  (3.6)  with  the  use  of  a  rotation  matrix  : 


r 

- 

h(a,<M,«)  *  a 

cos 0  -sin^ 

cost 

sin0  cor 4 

jsint 

• 

(3.7) 


The  two  descriptions  of  the  polarization  vector  h  of 
(3.3)  and  (3.7)  are  equivalent,  except  that  in  the  case  of 
(3.7),  the  antenna  is  descriped  by  geometric  parameters.  The 
parameters  a,  ^  t,  a  are  related  to  a^,  a^,  6^  and 


according 


oo  the  relations  (Appendix  A) 


a2-’a  *  +  a  1 

X  y 


2a  a 
x  y 


a  2  -a 

x  y 


tan2d 


2 


cor  5 


(3.8) 


PAGE  19 


sin2x  =  -  sinfi 

a  2+a  2 
x  y 

3.2.3  The  Stokes  Parameters 

A  practical  way  of  representing  the  state  of  polarization 

of  an  elliptical  wave  is  through  a  set  of  parameters,  which 

all  have  the  same  physical  dimensions.  Such  are  the  Stokes 

parameters  and  they  were  first  introducted  by  Stokes  in  his 

studies  of  partially  polarized  light.  These  parameters 

constitute  a  column  vector  g  (Stokes  vector)  and  they  can  be 

defined  in  terms  of  the  electric  field  vector  components,  or 

the  set(a  ,  a  ,  6  ,  5  )  or  (a,  0,  x,  a)  as  defined 
x  y  x  y 

in  Sections  3.2.1  and  3.2.2,  in  the  following  manner 
(Appendix  B)  : 

S0  =  lhJ2+lhyl2  =  Qx2  +  ayl  =  a2  =  I 

g.  =  | h  |  2“ |h,  |  2  =  a  2  -  a  2=  a2cos2,xcos2< $  -  Q 
l  x  y  x  y 

g„  *  2Re{h  h  }  =  2a  a  cosfi  =  a2cos2xsin2$  =  U  (3.9) 

*  x  y  x  y 

g„  =  -2Im{h  h  }  -  2a  a  sinfi  -  a2sin2.x  =  V 

3  x  y  x  y 

where  : 

gQ2  =  g12+S22+832  =  I2  =  Q1!+U2+V2  (3.10) 

The  component  g^  describes  the  intensity  while  g^,  g2 

and  g^  represent  the  polarization  of  an  ep  wave  since  they 

depend  on  the  orientation  of  its  ellipse  and  through  the 
ellipticity  angle  x  on  the  sense  in  which  the  ellipse  is 
being  described.  The  four  Stokes  parameters  have  the 
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dimensions  of  intens ity (power) ;  each  corresponds  not  to  an 
instantaneous  intensity  but  to  a  time  -  averaged  intensity, 
where  the  average  is  extended  over  an  interval  long  enough 
to  permit  practical  measurement.  Although  its  components  are 
physically  real  parameters,  the  Stokes  vector  is  a 
mathematical  vector,  i.e.  it  is  not  defined  in  a  three  - 
dimensional  physical  space  but  in  a  four-dimensional 
mathemathical  space.  The  description  of  polarization  through 
the  Stokes  parameters  is  widely  applicable  since  it  covers 
the  completely,  partially  or  unpolarized  waves.  The 
physical  proporty  of  the  4-corapcnent  vector  will  be 
discussed  in  another  forthcoming  report  dealing  with  quasi- 
coherent  pan-chromatic  wave  interaction. 

In  practice  the  modified  Stokes  vector  g  is  often 

**  m 

used  with  components 

8m0  =  Id+Q)  -  IV  -  -  Ix 

=i(I-Q)  =  |hy|‘  *  V  -  :y 

and 

gm2=g2’  gm3=g3  aS  given  in  Eq>  (3l9) 

3.2.4  The  Poincare  Sphere  : 


(3.11) 


The  Poincare  sphere  concept  consists  of  mapping 
polarization  states  on  points  of  a  sphere  for  completely 
polarized  ep  waves.  It  constitutes  a  useful  vay  of 
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representing  polarized  ep  waves  as  it  will  be  shown  as  the 

present  treatment  unfolds.  The  polarization  state  of  an  ep 

wave,  which  either  is  described  by  the  polarization  vector  h 

or  the  Stokas  vector  g,  can  be  represented  on  the  Poincare 

sphere.  Thus  one  can  map  the  polarization  vector  h  with 

complex  components  h  and  h  ,  given  by  (3.3),  on  the 

x  y 

Poincare  sphere  by  using  the  auxiliary  complex  parameter 


h  -jh  1-jP 

x  J  y  J 

u  =  -  = -  (3.12) 

hx+jhy  1+jP 

where  P  is  the  complex  polarization  ratio  and  it  is  equal  to 


P  =  -  (3.13) 

h 

x 

Using  (3.12)  and  (3.13),  the  polarization  vector  h  can  be 
represented  uniquely  on  a  point  on  the  Poincare  sphere  with 
spherical  coordinates  (r,  0,  <p')  which  are  given  by 

(Appendix  C)  : 

r=|  h  |2  =  a2 

|u|2-l 

co-latitude  9  =  arccos{ - }  =  tt/2-2t  (3.14) 

|u|2+l 


Im{u) 

longitude  01  =  2<t>  -  -phase{u}  =  -arctan{ - ■} 

Re{u) 

It  should  be  noted  that  absolute  phases  cannot  be 
represented  on  the  Poincare  sphere.  In  case  the  polarization 
state  is  described  by  the  Stokes  vector  g,  it  is  easily  seen 
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from  Eqs.(3.9)  that  its  g^,  g^,  components  can  be 

regarded  as  the  cartesian  coordinates  of  a  point  P  on  a 
sphere  of  radius  gg,  such  that  (tt/2 -2x )  and  20  are  the 

spherical  angular  coordinates  of  this  point (Fig. 3 . 3) . 

According  to  the  above  properties,  the  various  states  of 
polarization  can  be  mapped  on  the  Poincare  sphere  which  has 
the  following  properties (Fig. 3 .4)  : 

1.  The  Equatorial (x-y)  plane  divides  the  sphere  into  the 
left-sensed  upper  hemisphere  where  t  is  positive  (0  <  T  5 
ir/4)  and  the  right-sensed  lower  hemisphere  where  t  is 
negative  (-ir/4  £  t  <  0). 

2.  All  linear  polarizations  (t  =  0°)  are  represented  on  the 
Equator  with  horizontal  polarization(H)  at  zero  longitude  0' 
«  20  *  0  and  the  vertical  polarization(V)  at  the  antipodal 
location  0*  =  20  *  ir 

3.  Left-handed  circular  polarization(LC,T=iT/4)  is  mapped  on 
the  zenith  (gj^g^Q and  right-handed  circular 

polarization  (RC.ts-tt/A)  on  the  nadir 

Cg^~g2~®  »  ®3~'8q ) • 

4.  Any  two  orthogonal  polarizations  h  (a,0,t,o)  and 
h  _1_(a_L>0+ir/2, -t  ,aj_)  are  mapped  on  antipodal  pcdnts. 

5.  Statistical  polarization  can  also  be  mapped  on  the 
Poincare  sphere,  thermal  radiation  produces  random 
polarization  states  uniformly  distributed  over  the 
polarization  sphere. 


3.2.5  Polarization  Charts  : 


Polarization  charts  provide  a  useful  tool  for  a 
simplified  representation,  of  the  state  of  any  polarization 
vector  on  a  two-dimensional  chart  instead  of  using  the 
three-dimensional  Poincare  sphere,  which  though  is  still 
preferable  in  case  the  XPOL  and  COPOL  null  characteristics 
need  to  be  studied,  in  detail.  Deschamps (1951 , 1953)  in  his 
tutorial  study  on  tne  Hyperbolic  Protractor  and  Rumsey(1951) 
showed  how  various  mappings  on  the  sphere  can  be  related  to 
the  power  impedance (Smith)  chart  which  was  further  followed 
up  by  Huynen(1960) ,  Poelman(1971) ,  and  also  in  the  Russian 
literature  as  e.g.  in  Kanareykin  et  al(1966 , 1968) .  There 
exist  many  types  of  polarization  charts  and  here  we  will 
briefly  review  the  properties  of  some  charts  which  are  being 
used  most  frequently. 

1.  Rumsey  p-and  q-Charts  : 

Rumsey(1951)  has  used  the  impedance  concept  which  reduces 
the  field  problem  to  a  transmission  line  problem.  The 
analysis  of  impedance  transformation  or.curing  in 
transmission  line  theory  can  again  be  simplified  by  working 
in  terms  of  reflection  coefficients.  The  impedance  concept 
is  successful  in  such  applications  largely  because  the 
tangential  components  of  the  electric  and  magnetic  fields 
are  continuous  at  a  surface  separating  two  different  media. 


For  the  same  reason,  he  used  the  polarization  ratio  P  defind 
as  the  ratio  of  two  orthogonal  tangential  components  of 
electric  field,  to  be  equally  valuable  in  the  analysis  of 
polarization. The  polarization  ratio  P  is  analogous  to 
impedance,  then  Rumsey  finds  that  the  analogue  of  the 
reflection  coefficient  is  the  ratio  of  the  right-handed  and 
left-handed  circularly  polarized  components  which  are 
equivalent  to  the  linearly  polarized  components  used  to 
define  P.  Following  Rumsey(1951)  : 


circularly  polarized  components.  Note  that  p=jP 
=j (Polarization  ratio).  The  transformation  from  p  to  q  is 
thus  identical  to  the  transformation  from  the  current 


reflection  coefficient  to  the  normalized  impedance.  In  view 
of  the  symmetry  of  the  transformation  we  can  think  of  p,  and 
vice  versa  q,  as  reflection  coefficient  or  impedance.  By 
using  these  relations,  Rumsey  has  developed  the  q  and  p 
charts  for  representing  any  polarization  vector. 


Fig. 3. 5  shows  how  the  orientation  and  shape  of  the 
polarization  ellipse  are  represented  using  th.  Smith- 
Buschbeck  Chart  as  the  q-plane.  Fig. 3. 6  illustrates  the 
representation  using  tue  Carter-Schmidt  diagram  as  the  p- 
plane. 


Fig.  3.5  :  q-plane  Polarization  Chart  (Rumsey  1951). 
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For  points  on  the  q-plane  Fig. 3. 5,  (Rumsey  (1951))  : 

a)  We  have  right-handed  polarization  for  all  points  irside 
the  unit  circle,  and  left-handed  polarization  for  all  points 
outside. 

b)  The  origin  represents  right-circular  polarization,  and 
the  point  at  infinity  represents  left-circular  polarization. 

c)  Constant  axial-ratio  contours  are  identical  to  the 
circles  of  constant  dWR  on  a  Smith  chart.  The  axial  ratio  is 
equal  to  the  SWx  obtained  by  treating  q  as  a  reflection 
coefficient . 

d)  A  point  on  the  unit  circle  represents  linear  polarization 
at  an  angle  equal  to  one-half  the  polar  angle. 

a)  The  angle  <f>  between  the  x  axis  and  the  major  axis  of  the 
polarization  ellipse  is  related  to  the  polar  angle  0  on  the 
plane  by  the  simple  relation  2^=0. 

f)  the  locus  of  the  points  representing  polarizations  for 
which  Ih^/hyl  is  constant,  is  the  set  of  circles  passing 

through  the  points  q=±l  (the  short  and  open-circuit  points), 

the  orthogonal  set  of  circles  is  obtained  if  the  phase  of 

h  /h  is  constant, 
x'  y 

For  points  on  the  p-plane,  Fig. 3. 6,  (Rumsey  (1951))  : 
a)  Points  in  the  right  half-plane  represent  right-handed 
polarization  and  points  in  the  left-plane  represent  left- 
handed  polarization. 
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b)  Points  on  the  vertical  axis  represent  linear  polarization 
at  various  angles . 

c)  Circular  polarization  is  represented  by  the  points (1,0) 
and  (-1,C;. 

d)  Constant  axial-ratio  contours  are  identical  to  the 
circles  of  constant  standing  wave  ratio (SWR)  on  the  Carter- 
Schmidt  impedance  chart.  The  axial  ratio  is  equal  to  the  SVHR 
obtained  by  treating  p  as  an  impedance. 


Fig.  3.6  :  p-plane  Polarization  Chart  (Rumsey  1951). 
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2.  Deschamps  0  and  S  charts  : 

Deschamps(1951 , 19S3)  has  used  the  polarization  ratio  P 
which  can  be  deduced  from  the  parameters  <A,  x  (orientation 
and  ellipticity  angles),  to  derive  the  two  sets  of  equations 
cos2t  *  cos2rcos20 
tan  6  »  tan2tc»c2$ 
or 

tan2^  =  tanSYcosS 
sin2t  *  sin2y«infi 

where  P  «tanlTe^  (5  is  the  phase  difference  between  y  and 
x  channels  and  tanJT  is  the  amplitude  of  the  polarization 
ratio).  He  has  used  two  projections  on  the  equator  of  the 
sphere.  The  first  one,  an  orthographic  projection  (chart  0) 
is  shown  in  Fig.  3. 7.  Tire  other,  a  stereographic  projection 
from  the  nadir,  is  shown  on  Fig. 3.8  (Chart  S). 

In  both  cases  by  drawing  the  lines  along  which  ^  and  t 
are  constant  (meridians  and  parallels)  and  the  lines  along 
which  6  and  V ,  or  the  ratio  tanJf,  are  constant,  we  have  a 
method  of  direct  conversion  from  any  two  of  these  parameters 
to  the  others . 


The  lines  5  constant  are  circles  on  Chart  S  and  ellipses 
on  Chart  0.  The  lines  ?'  constant  are  straight  lines 
perpendicular  to  HV  on  Chart  0  and  circles  on  Chart  S. 


ituiibiL 
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3.  Huynen's  Chart  ; 

Huynen(1970)  has  presented  two  polarization  charts,  each 
representing  one  hemisphere  (top  and  bottom) ,  to  map  the 
whole  sphere  on  a  plane.  Fig. 3. 9  shown  such  a  circular 
polarization  chart,  which  maps  all  positive  or  left-sensed 
polarizations.  The  circumference  of  the  circular  chart  gives 
all  linear  polarizations  and  the  chart  is  left  circular. 
Notice  the  effect  of  2$  on  points  of  the  chart  such  that 
"horizontal"  polarization  ($“Q°)  is  mapped  on  the  extreme 
right-  hand  side  of  the  chart,  while  "vertical"  polarization 
(2^“180°)  is  mapped  on  the  extreme  left-hand  side.  All 
points  on  the  chart  represent  polarization  with  orientation 
$“45°.  Note  the  intereating  fact  that  the  radial  distanco  of 
a  point  on  the  polarization  chart  is  measured  by  cos2t  . 


X  , . -<K 

/  X  T  > 

/  'Q)  V?v"S 


i - 

XIX  J  , 

"r"  X 

' - - - ^ 


Fig.  3.9  :  Huynen’s  Polarization  Churt  (Huynen  1970) 
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4.  Poelman'a  Modified  Chart  : 

Poelman(1971)  has  modified  Huynen's  chart  which  was 
explained  before.  He  used  in  his  modified  polarization  chart 
the  polar  coordinates  20  and  r  such  that  : 

r  =*  tant  *» 

as  then  the  geometric  parameters  0  and  r  can  easily  be  read 
off.  Note  that  Poelman  used  r=tant  instead  of  cos2t  which 
was  used  by  Huynen. 

In  Fig. 3. 10  the  distinguished  polarizations  with  the 
proposed  code  numbers  for  the  right-sensed  polarization  are 
presented  on  the  modified  polarization  chart.  A  similar 
chart  for  the  left-sensed  polarization  can  be  given,  where 
for  the  linear  polarization  the  same  code  numbers  can  be 
taken  as  given  in  Fig. 3. 10. 

In  Table  3.1,  the  149  used  polarizer  settings  (0,6)  for 
the  right-sensed  polarizations,  including  the  28  linear 
polarization  and  the  corresponding  code  numbers  are  given. 
It  follows  that  a  group  of  270(149+149-28)  different 
polarizations  (right-and-left-sensed)  is  chosen  to  represent 
all  posoible  radar  polarizations. 

In  Table  3.2  the  polarization  characteristics  (r  and  0) 
and  corresponding  code  numbers  are  listed  for  the  right- 
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Tab;3  3*1:  Survey  of  the  used  right-sensed  polarizations;  the  polarization 

characteristics  corresponding  with  the  1 6 1  code  numbers  are  presented 
in  Table  3-2,  and  the  areas  they  represent  on  the  polarization  chart 
are  giver  in  Fig.  3.10  (Poelman,  1971). 
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code 

No. 

H 

$ 

code 

No. 

r 

$ 

code 

No. 

r 

$ 

1 

D 

0° 

19 

0 

115.5° 

37 

0.11 

51.7° 

B 

6. 5° 

20 

0 

122° 

38 

0.  10 

58.3° 

3 

0 

21 

0 

128.5° 

39 

0.09 

64.9° 

4 

0 

19-5° 

22 

0 

135° 

40 

0.07 

70.9° 

5 

0 

25-5° 

23 

0 

141.5° 

41 

0.10 

78.2° 

6 

0 

32° 

24 

0 

148° 

42 

0.07 

84.9° 

7 

0 

38.5° 

25 

0 

154.5° 

43 

0.11 

90° 

8 

0 

45° 

26 

0 

160.5° 

44 

0.07 

95.1° 

9 

51-5° 

27 

0 

167° 

45 

0.10 

101.8° 

10 

0 

58° 

28 

0 

173.5° 

46 

0.07 

109. 1° 

11 

0 

64-5° 

29 

0.11 

0° 

47 

0.09 

115.10 

12 

H 

70.5° 

30 

0.07 

5.1° 

48 

0.10 

121.7° 

13 

B 

77° 

31 

0.10 

11.8° 

49 

0.11 

128.3° 

14 

0 

83*5° 

32 

0.07 

•19.1° 

50 

0.11 

135° 

15 

1  : 

90° 

33 

0.09 

25.10 

51 

0.11 

141.7° 

16 

96.5° 

34 

0,10 

31.7° 

52 

0.10 

148.3° 

17 

103° 

35 

0.11 

38.3° 

53 

0.09 

154.9° 

10 

0 

109.5° 

36 

0.11 

45° 

54 

0.07 

160.9° 

Table  3,2:  Survey  of  the  polarization  characteristics  ellip- 
ticity  ratio,  r,  and  orientation  angle,  cor— 
responding  with  the  code  numbers  for  the  right- 
sensed  polarizations.  (Poelman,  1971) 
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Fig.  3.11 
1977,73). 
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Poelman’s  modified  Polarization  Chart  (b)  (Poelman 
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polarization,  including  the  linear  polarization. 

Poelman£ 1977, 1978)  has  used  also  another  version  of  the 
modified  polarization  chart  Fig. 3. 11  which  is  very  similar 
to  the  one  shown  in  Fig. 3. 11  but  with  using  the  polar 
coordinates  20  and  p=l-r. 


5.  Optimal  Polarization  Charts  : 

The  optimal  polarization  charts  (Kennaugh,  1952,  Huynen, 

1970)  will  be  discussed  in  the  forthcoming  reports.  On  these 
charts,  the  C0P0L  and  XPOL  nulls  will  be  drawn  for  different 
aspect  angles  and  different  frequencies  for  various  targets 
and  sea  clutter. 

3.2.6  Canonical  Polarization  Pairs 

The  power  PCh^.lu)  received,  by  a  receiving  antenna 

with  polarization  -j (aj >®j > $ j »Tj )  due  to  that 

transmitted  by  h . (a . ,a . ,0 . ,T , )  is  found  to  be 
“1  1  11  1 

given  in  general  by  (Huynen  1970  ): 

P(h.,h.)  =  ^a . 2a . 2 [ l+sin2x ,sin2t ,+cos2 (0 ,+0 . )cos2i . cos2r , ]  (3.15) 

-j  i'  *  l  j  1  j  i  j  .i  j  iJ 

Thus  by  selecting  the  parameters  of  the  polarization  pairs 

(hj,h^),  one  can  determine  according  to  (3.15)  the 

magnitude  of  the  received  power.  There  are  four  such 
canonical  polarization  pairs  (Fig. 3. 12)  of  special 
importance  since  they  give  rise  to  situations  frequently 
occuring  in  practical  transimission  reception  interactions. 

These  are  : 


ttiiMifeitiiittii 
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ORTHOGONAL  :  h.  =  {  h.}*, 

-  -j  — i 

h  ,=h.(a.  and  a.  arb.  .=$  ,+iif.T  =-t , 

-orthogonal  -jj  j  j  i  j  i 

antipodal ,  but  undetermined  in  a  and  a: 

CANNOT  BE  USED  AS  A  UNIQUE  DESCRIPTOR. 


TRANSVERSE  :  h ,=h . (-a . . -d  . )  :  PCh^h.^P  =a.*a,1 
-  -j  -iv  i»  Ti'  j  * — i-'  max  i  j 


h  =h ,  (a  ,=a.  .«  .--a ,  =-$  ,  ,t  =t  , ) 

-transverse  -j  j  i’  j  i,rj  i*  j  i' 

OPTIMAL  RECEPTION  :  ANTENNA  MATCHING. 


SYMMETRIC  :  h  (=h,  (-*  ,  -t  ,  )  :P(h .  ,h .  )=a , **cos22t  . 

-  -j  -i  Ti  i '  t  *— i  i  i 

h  .  .  =h  ,  (a  ,=a,  ,a  .=<x,  ,x  =-t .  ) 

-symmetric  -j  j  i’  j  i’  j  i‘j  r i' 

MOST  FREQUENCY  TARGET  POLARIZATION. 


CONJUGATE  : 


h  .  .  =h,  (a=a,  ,o  .=*-a,  .  ,x  =-t  , ) : 

-conjugate  -j  j  i *  j  i’  j  i’  j  i' 

PRECIPITATION  (CIRCULAR  CLUTTER')  . 


orthogonal 


SYMMETRIC 


TRANSVERSE 


CONJUGATE 


Fig.  3.12  :  Basic  Polarization  Pairs. 
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3.3  Scattering  MatrixfSI  ; 

The  general  radar  target  detection  technique  consists  of 
having  an  electromagnetic  radiation  originating  from  a  radar 
transmitter  scattered  by  the  object  under  detection  and 
sampled  by  the  radar  receiver.  The  measure  of  the  intensity 
of  the  scattered  radiation  in  the  far  zone  from  the 
scattering  obstacle  was  usually  described  by  the  scalar 
radar  range  equation  expressed  in  terms  of  the  radar  cross- 
section  (RCS)  a  given  by  : 

!Es(es,*g)l2  |Hs(es,*s)l2 

<J=lim4irR* - 3lim4irRa -  (3.16) 

R"“  lEi(9i^i)|*  |Hi(ei,*.)|a 

where  R  is  the  range  (distance  from  target  where  the 

scattered  radiation  is  observed,  Fig. 3. 13),  E  S,  H  s  are 

the  scattered  electric  and  magnetic  field  vectors  at  the 

observation  point  at  direction  (0  ),  and  similarly 

s  s 

e\  define  the  incident  electric  and  magnetic  field 
vectors  with  direction  angles  0^,  0^. 


SOURCE 


RECEIVER 


Fig.  3.13  :  Illustration  of  Radar  Range  Equation. 
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Eq.(3.16)  gave  unsatisfactory  results  with  clutter 
precipitation  present  thus  yielding  limited  radar  target 
detection.  Recognizing  these  difficulties,  Sinclair (1948) 
showed  that  by  introducing  polarization  in  the  radar  range 
equation  the  problems  could  be  overcome.  He  proved  that  a 
radar  target  acts  a  polarization  transformer  and  this 
polarization  transformation  can  be  described  by  a  matrix 
known  as  the  scattering  matrix.  We  note  that  in  Optics 
Jones(1941)  introduced  a  similar  matrix(Shurclifii,  1962). 


s  i 

Let  the  polarization  vectors  h  and  h  for  the 
scattered  and  incident  radiation,  respectively,  be  defined 
in  terms  of  general  orthogonal  polarization  base  vectors 

A  A 

» hg )  i *  e • 


V 

and 

h1* 

v. 

t  i 

V. 

The  complex  components  h^,  h^  can  represent  any 


polarization  by  relative  magnitudes  and  phases.  The 

scattering  matrix  transforms  the  transmitted  polarization 

i  s 

h  into  the  polarization  of  the  scattered  field  h 

according  to  i 


hS(es#*s)=t  m9;J,*s,9i,ti)]hLWi,*i) 

where  [  /a]  in  case  *of  A=H,  B=V  is  the  "linear  polarization" 
restricted  scattering  matrix  with  -absolute  phase,  i.e. 


^lOgge-^BB'W 


(3.17) 


f J. 18) 


(3.19) 


The  scattering  matrix  represents  the  radar  target  for  a 
given  frequency  and  fixed  aspect  angle  .  Thus  at  these  fixed 
frequency  and  aspect  angle,  the  target  is  described 
completely  by  the  scattering  matrix.  However,  in  special 
cases  one  and  the  same  matrix  may  represent  properties  of  a 
set  of  different  targets  (Huynen,  1970).  A  particularly 
useful  feature  of  the  scattering  matrix  representation  is 
that  its  intrinsic  properties  are  a  function  of  the  target 
configuration  and  do  not  depend  on  the  measurement  technique 
or  measuring  equipment.  However,  the  scattering  matrix  with 
absolute  phase  does  depend  on  the  target  displacement  along 
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ml  ini'. .  .  ,!"■ 


the  line  of  sight.  In  order  for  the  *■  ® J SMA  matri>t:  t0 

specified  completely,  eight  reel  numbers (four  magnitudes  and 
four  phases)  have  to  be  measured.  The  measurement  of  the 
absolute  phase  of  a  target  becomes  very  formidable  because 
of  its  dependence  on  the  location  of  the  target,  the 
direction  of  illumination,  surface  configuration  and  the 
radar  frequency.  Hence,  one  has  to  use  highly  elaborate 
measuring  techniques.  Moreover,  the  absolute  phases  cannot 
bo  mapped  on  the  Poincare  sphere  (Thiel  1970)  and  thus  the 
^SMA  raatr*x  cannot  have  a  unique  representation  on  the 

sphere.  Beccuse  of  these  difficulties,  the  scattaring  matrix 
with  relative  phase (SMR),  whore  any  one  of  the  phases  of  the 
matrix  elements  can  be  set  to  an  arbitrary  constant  (moat 
commonly  the  phase  of  the  cross -polarized  term  i.e. 

is  chosen  to  be  zero) .  Thus  the  scattering  matrix  with 
relative  phase  [S)g^  is  according  to  (3.22)  ; 


i 


*  ‘SMR 


SAA  SAB 


SBA  SBB 


with  *AA'*BA‘W0>  Und 


*AB“° 


(3,23) 


This  simplification  reduces  the  number  required  to  specify 
completely  the  scattering  matrix  from  'Sight  to  seven  real 
numbers . 

Another  important  feature  of  defining  relative  phases 
instead  of  the  absolute  ones  is  that  relative  phase  terras 


III.  ,li,.,..’tllnj  lliJnl;,Mm. Ill,,,1  iJn  <:i,..l..  II  ■illiJi  iM  „:il'  il|l| fulfil  illl.lJik'lM'/LU.'iki^'Mlik.jljl.-j.f  i1illi:l',i,'fl,'j ,  lilJli'ilJL^'l'  .  II  ,,l  .■  .  ......  .  :  Lj y '.f ..if  ||, £ 
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can  be  recovered  from  amplitude  only  data  (Kennaugh, 1949) . 
la  the  case  of  "blstatic"  scattering  one  has  :  8  ^8,, 

S  X 

^AB^BA*  *n  "monosta,tic"  caae  (both 

transmitter  and  receiver  antennas  are  at  the  same  location), 
we  have  because  the  reciprocity  theorem  and 

conservation  of  energy  should  be  satisfied  for  propagation 
in  an  isotropic  medium,.  We  note  that  for  anisotropic 
scatters  this  condition  need  not  be  satified,  Thus  the 
scattering  matrix[S becomes  symmetric  in  the 

"monos t at ic"  case  and  only  five  real  numbers  (three 
amplitudes  and  two  phases)  are  needed  to  specify  it 
completely. 


Using  (3.17),  (3.21)  and  (3.22),  we  have 


*•*  M 

*• 

r  jl 

V 

m 

SAA 

SAB 

V 

SBA 

SBB 

M»  Ml 

• 

w 

lUioJlLkt'.ii'.lL  *,y  •/ii-.id.lliL. :AL 


The  polarisation  of  the  electromagnetic  wave  in  (3.24)  is 
defined  with  respect  to  two  orthogonal  components  which  are 
in  general  elliptical.  The  two  special  cases  of  elliptical 
polarization  i.e.  linear  and  circular  polarizations  are  the 
most  commonly  used  for  reference  systems.  If  we  represent 
the  linear  polarization  vector  with  (hx>hy)»  and  the 

circular  polarization  vector  with  (h  ,h.)  as  its  light 

r  JC 

and  left:  circular  components,  respectively,  then  the 
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relations  between  the  scattered  fields  in  the  two 
polarizations  are  given  by  (Long  1966)  : 


\  =  i(h+jhy),  hr  =  i(hx-Jhy) 


(3.25) 


The  scattering  matrices  for  the  two  polarizations  are 
related  through  the  formulas  : 


n  =  |4(S  -S  )+JS  |  ,  S  =S 
It  1  xx  yy'  J  xy1  ’  xy  yx 

c,  -  |i(s  +s  ) I ,  c.  =c  „ 

£r  1  xx  yy  1 ’  Hr  r l 

C  =  |i(S  -S  )-jS  | 
rr  1  xx  yy  xy1 

where  the  are  the  elements  of  the  scattering  matrix 

with  respect  to  the 

circular  polarization  basis  and  are  the  corresponding 
elements  with  respect  to  linear  polarization  basis. 


3.4  Scattering  Matrix  Transformation  Invariants  : 

3.4.1  The  unitary  transformation  matrix  [T1  ; 

Assuming  reciprocity  holds,  there  exists  an  infinite 
number  of  general  pairs  of  orthogonal  elliptical 

/N 

polarization  basis  vectors  h^.hg  and  an  infinite  number 

of  possible  invariant  transformations.  Numerically,  the 
transformation  properties  of  the  scattering  matrix  [S(A,B)] 
(assuming  no  polarization  losses)  from  one  orthogonal  pair  : 


(3.26) 


-  ~  hAhA  +  hBhB 


(3.27) 


to  another  orthogonal  pair  : 


h  =  hA*v  ♦  ^ 


(3.28) 
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can  be  transformed  to  another  new  scattering  matrix 
[ S ' (A ’ , B * ) ]  by  using  a  unitary  transformation  matrix  [T]  . 
This  matrix  [T]  relates  the  polarization  vector  components 
in  the  two  orthogonal  basis  sets  (A,B)  and  (A',B').  To  find 
the  matrix  [T] ,  the  orthogonality  relationship  between  the 
components  of  each  of  the  two  orthogonal  basis  vectors 
should  be  used  e.g. 

A  A  *jV 


h  *h.  =  fi. . 

i  J  ij 


where 


5ij 


is  the 


Kronecker 


delta 


f unct ion 


(3.29) 


(6^=1, i3j, 6^=0, i/j)  and  (i,j)  may  be  equal  to  A  or  B 
in  the  first  basis  set  and  A*  or  B1  in  the  second  set. 


From  (3.27)  to  (3.29),  we  have  : 


A  *  .  *  .A,** 

-*hA  =  hA  *  hA  (hA  -hA  )  +  WhA  > 


(3.20) 


and 

-*hB  hB  “  hA  ^A  *hB  ^  +  hB  ^rtB  *^B  ^ 

By  rewriting  (3.30)  and  (3.31)  in  matrix  form,  one  has 


A 

N 


”  .  A  Vf  A  A  * 

hA  *hA  hB  *hA 


A  ,  A.  *  A  ,  A  * 

LhA*hB  hB*hB. 


h.  ’ 
A 


V 


(3.31) 


(3.32) 


or. 


h(A,B)  =  [T]  h(A'.B') 
where  : 


■A  A  *  A  A  * 

hA  *hA  hB  *hA 


(3.33) 


[T]  = 


(3.34) 
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To  determine  the  matrix  [T] ,  the  relationship  between  the 
two  sets  of  the  orthogonal  basis  vectors  should  be  known. 

For  example,  let  : 

hA'  =  otjh^  +  a  2^8  (3.35) 

and 

V  =  +  &2*B  (3-36) 

These  two  new  unit  vectors  h^' ,  h^'  in  the  orthogonal 

basis  (A1, B')  should  satisfy  the  orthogonal  1 ty  relationship 
of  (3.29).  For  (3.29)  to  hold,  the  following  must  be  met 


laj1  +  |«2|a  *  1,  (3.37a) 

IBJ2  +  |  32!  2  «  1,  and  (3.37b) 

°1B1*  +  *2*2  =  °  (3.37c) 

By  substituting  (3.35)  and  (3.36)  into  (3.34)  and  using 
(3.29),  the  matrix  [T]  can  be  written  in  the  form  : 


(3.38) 


This  means  that  there  exists  an  infinite  number  of 
unitary  matrices  [T]  given  by  (3.38),  which  satisfy  the 
conditions  (3.37).  For  example,  if  we  choose  : 

^1  =  '“2**  &2  =  “l*  (3.39) 

which  satisfies  (3.37b  &  c)  provided  that  (3.37a)  is 

satisfied,  Eq.(3.38)  yields  : 
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h^A.B)  =  [T]  h1  C A '  ,B '  )  (3.44) 

and 

hS(A,B)  =  [T*]hS(A' ,B')  (3.45) 

* 

We  use  [T  ]  instead  of  [T]  in  (3.45)  in  order  to  preserve 
the  same  sense  of  polarization  with  respect  to  both 
coordinate  systems  of  the  incident  and  scattered  radiation 
(Graves  1956,  Maffett  1968).  From  (3.45)  : 

hS(A',B')  =  [T*  ]  ~  1h_S  (A ,  B )  =  [TT]hS(A,B)  (3.46) 

Using  (3.21)  and  (3.44)  one  has  : 

hS(A,B)  «  [S(A,B) ]hi(A,B) 

*  iS(A,B)][T1hi(A,3B')  (3.47) 

Substituting  (3.47)  into  (3.46),  we  obtain  t 

hS(A',B')  *  [TT][3(A,B)HT]hi(A\B’) 

=  [S,(A,,B,)]hi(A',B')  (3.48) 

where  : 

t S  * (A  * , B  * ) ]  =  [TT][S(A,B)][T]  (3.49) 

From  (3.49)  and  (3.42),  we  obtain  : 

S'A,A,  =  (l+pp*)‘1[SAA+p*SBB+p(SAB+SBA)l 

S'A,B,  =  a+PP*)'li-P*SAA+pSBB+SAB-pp*SBA]  (3.50) 

S’b'A'  =  (l+PP*)'li-P*SAA+PSaB+EBA'pp*SABi 

s'B,B,  =  a+pp*)'1tp*!sAA+sBB-p*csAB+sBA)] 

Using  (3.50)  we  find  that  the  det{ [S(A,B) ] }  and  the 
span  {[S(A,S)]}  are  transformation  invariants  (Appendix  D) , 
according  to  : 
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det{ [S(A,B) ] }  =  det{ [S ' (A* ,B ' ) ] }  =  invariant,  (3.51) 

and 

span{  [S(A,B)  ]  }  =  I  SMI2+I  SABI  J+|  SBA|2+|SBB|a  =p 
=  span{ [S ' (A* ,B' ) ] ) 

=  !sVA*l2+iSVB-l2+|S,B*A'l2+IS,B*B*la 

3  invariant  (3.52) 

We  note  that,  if  SAfi  =  S^,  then  S ' A , ,  =  S ' , A ,  , 


i.e.  ,  if  reciprocity  is  satisfied  for  any  one  pair  of 
orthogonal  polarizations,  it  is  satisfied  for  all  such 
pairs.  Furthermore,  we  must  emphasize  the  important  property 
that  for  any  one  given  aspect  and  for  one  frequency,  the 
transformation  occurs  on  one  and  the  same  polarization 
sphere  of  radius  p  3  span{ [S(A,B) ] }  3  span  { [S ' (A' ,B ' ) ] ) . 
Thus, if  (S(A,B)]  is  known,  the  new  scattering  matrix 
£ S * (A * , B * ) ]  for  any  other  orthogonal  base  (A',B')  can  be 
determined  if  the  relationship  between  the  two  sets  of  bases 
(A,B)  and  (A',B')  are  known.  This  is  shown  for  example  in 
the  transformation  from  linear  to  circular  polarization 
basis  vectors  in  (Long  1966).  In  case  of  polarization 
losses,  the  properties  of  the  coherency  matrix  need  to  be 
used  (Kraus  1966),  and  transformation  will  not  occur  on  the 
same  polarization  sphere  (Deschamps  1951),  as  is  discussed 
in  (Thiel  1970)  and  will  be  further  analyzed  in  one  of  our 
forthcoming  reports . 
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3.5  The  Optimal  Polarizations  : 

It  was  first  shown  by  (Kennaugh  1949-1954),  that  there 
exist  two  pairs  of  optimal  polarization  which  can  be 
associated  with  (3.50)  and  are  useful  to  express  the  five 
independent  real  components  of  the  scattering  matrix  [S]  in 
the  monostatic  relative  phase  case  on  the  polarization 
sphere  of  radius  p  subject  to  (3.52). 


The  CO -POLARIZATION  (COPOL)  NULL  PAIRS  for  minimal 
polarization  are  obtained  from  (3.50)  by  setting  S ' ^ A * 

and/or  S  * „ , _  t  to  zero  so  that  in  the  bistatic  case  : 

D  D 

‘‘W*  ^SAB+SBA>‘-4SAASBB 


CO 


1,2 


2S 


BB 


reducing  for  the  monos*- atic  case  to  : 


■SAB*t/SAB2-SAASBB 


CO 


BB 


The  CROSS -POLARIZATION  (XPOL)  NULL  PAIR  for  maximal 
polarization  is  obtained  from  (3.50)  by  setting  S'^,gi  or 

3 '  ,  .  .  to  zero  so  that  in  tue  bistatic  case  : 

B  A 


-b±  ^  b* -4ac 


1,2 


2a 


where  : 


a  "  SBBSBA  +SAA  SBA 


b  =  "^'SB2*  I‘SABaBA  +SAB  SBA''SAa'^ 


C  =  '(SBB  SAB+SAASAB  -1 


(3.53) 


(3.54) 


(3.55a) 


.55b) 

■  i*  j iWtv1  .". 


PAGE  S3 


reducing  for  the  monostatic  case  to  : 

b  =  "{iSBBl2"|SAAi1}  311(1  C  =  “a*  (3.55c) 

If  we  let,  in  general  : 

1-jp 

u  =  -  (3.56) 

1+jp 

with  p  being  the  transformation  parameter  defined  by 
Eq.(3.50),  the  coordinates  resulting  for  (3.53)  to  (3.55)  on 
the  Poincare  sphere  are  given  by 

|u|2-l 

Colatitude  :  8  =  arccos  -  (3.57) 

|u|a+l 

Im{u) 

Longitude  :  i $'  =  -arctan  -  =  -phase{u>  (3.58) 

Re{u} 

On  the  basis  of  these  expressions,  examples  of  calculations 
of  the  optimal  polarization  pairs  for  radar  targets  and 
clutter  are  given  in  Chapter  Four.  It  should  be  noted,  as 
is  illustrated  in  Fig.  (3. 14),  that,  for  the  monostatic  case 
the  XPOL  and  COPOL  nulls  lie  on  one  main  circle,  *'  it.  the 
XPOL  nulls  are  antipodal,  and  the  connecting  line  bisects 
the  great  circle  arc  between  the  COPOL  nulls.  This  means 
that,  if  we  determine  the  COPOL  nulls,  it  is  easy  to  find 
the  location  of  the  XPOL  nulls  but  not  vice  versa.  Also,  if 
one  XPOL  null  and  one  COPOL  null  are  given  (e. g.  by 
measurements),  the  location  of  the  other  two  nulls  can  be 
found.  By  using  the  COPOL  null  polarizations,  the  scattering 
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matrix  can  be  graphically  represented  by  two  points  on  the 
surface  of  the  Poincare  sphere  of  radius  p.  The  coordinates 
of  the  scattering  matrix  or  of  the  target  signature  and  the 
radius  of  the  sphere  represent  the  function  of  the  target 
reflectivity,  and  in  principle  contain  the  same  information 
as  [S]  for  the  monostatic  relative  phase  case. 

Details  of  the  intrinsic  properties  of  the  optimal 
polarization  nulls  on  the  Poincare  sphere  will  be  discussed 
in  one  of  our  forthcoming  reports  together  with  the  analysis 
of  model  generated  and  measurement  data. 

3.6  The  Stokes  or  Mueller  Matrix  and  its  Relationship  with 
the  Scattering  Matrix  : 

Whereas,  in  the  coherent  case  the  elements  of  the  2x2 
Sinclair  matrix  [S(A,B)]  are  additive,  in  the  incoherent 
case,  the  time-averaged  Stokes  parameters  of  the  4x4  Mueller 
matrix  are  additive.  The  Mueller  matrix  [M]  relates  the 
scattered  and  the  incident  Stokes  vectors  in  the  following 
way  : 

gS(A,B)  =  [MJg^A.B)  (3.59) 

where  the  Stokes  vector  g  =  (gQjg^g^gg)  = 

(I,Q,U,V)  .is  given  by  (3.9)  in  terms  of  the  phasors  of  its 
corresponding  polarization  vector  : 

±mbj?A%*B  (3-60) 

Similarly,  the  modified  Mueller  matrix  [Mm]  relates  the 
modified  scattered  and  incident  Stokes  vectors  through  the 


equation  : 
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anjS(A,B)  = 

where  !o  =  CiCI+Q)  ,i(I-Q)  ,TJ,V)  given  by  (3.11). 


(3.61) 


In  the  following  sections  we  will,  as  given  by  (3.11), 
derive  the  relationship  between  the  matrices  [M] ,  [Mm]  and 
[S].  In  the  first  subsection,  the  derivation  of  the  elements 
of  [M]  and  [Mm]  from  [S]  is  given,  while  in  the  second 
subsection,  the  inverse  problem,  i.e.  ,  the  derivation  of  the 
elements  of  [S]  from  [M]  or  [Mm]  is  solved.  Also  the 
relation  between  [M]  and  [Mm]  is  given  in  the  third 
subsection. 


3.6.1  Derivation  of  the  Mueller  [Ml  or  Modified  Mueller  [Mm] 
Matrices  from  the  Scattering  Matrix  fS]  : 

a)  Derivation  of  the  Mueller  Matrix  [Ml  from  the  scattering 
matrix  [S]  : 

The  scattering  matrix  [S(A,B)]  is  given  in  the  orthogonal 
basis  (A,B)  by  (3.22)  which  in  the  bistatic  case  has  four 
complex  quantities  or  it  possesses  eight  real  quantities  for 
the  absolute  phase  case.  The  complex  quantities  are  S^, 


AB’ 
S 


*BB‘ 


SBA  ***  SBB  ’Whil® 

'  *  'sab'*  ^sba'’  'sbb'’ 

In  the  monostatic  case. 


real  quantities  |S^|,  |S^gj 


*SB  because  iSBA!  =  ^AB1  md 


the  real  quantities  are 

*AA’  *AB’  *BA* 
we  have  =  S^g  or  six 

*  'sbb'*  *AA’ 

*'BA  =  *AB-  In  this 


section  the  elements  of  the  Mueller  Matrix  [M]  for  both  the 
bistatic  and  monostatic  cases  are  derived  from  the  elements 
of  the  scattering  matrix  [S(A,B)],  as  shown  in  Appendix  (E) 
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i)  Blstatic  Cass  : 
°11  =  * 
m12  =  * 

°13  =  R 

*14  =  1 
“21  =  * 

“22  =  * 


1m|,','|sba;,+|sabi 

2  +ISBB 

!*aI*+Isba'*-'sab 

1 2'Isbb 

'aasab  +sbasbb  } 

'AASAB  +SBASBB  } 

;  1 *- 1 s  |*+|s  I 

1  -Is 

AA 1  1  BA1  1  AB' 

1  BB 

!AAl2"iSBAl2‘iSABl 

'*  +ISBB 

m23  =  Re{SAASAB  "SBASBB  } 


m24  Im{SAASAB  “SBASBB  } 


m31  =  Re{SAASBA*+SABSBBA) 


m32  Re*SAASBA  “SABSBB  } 


“*■-  =  Re{SAASB3  +SABSBA  * 


33 


m34  =  Im{SAASBB  "SABSBA  } 


*41  =  "Im^SiASBA  +SARSRR  ^ 


AA  BA  AB  BB 


*42  =  "Im{SAASBA  "SABSBB  } 


~  '^^AA^na  +®ad^ija  ) 


AA  BB  AB  BA 


m44  =  Re(-SAASBB  ’SABSBA  } 


li)  Monostatic  Case  : 

*11  =  ^lSAA'2+2|SABi2+|SBB|l} 


*12  =  ^ISAa'2-ISBbI2} 


l2> 


l2> 

I2) 


m,,  =  Re(SAAS.n  +SAnSDD  > 
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°14  "  Im(SAASAB  +SABSBB  } 


m21  =  °12 


m 


22 


"23  =  Ee(SAASAB*-SABSBB*> 


“24  * 


ImfSAASAB  'SABSBB  1 


m31  =  m13 


m32  “  m23 


m33  = 


Ro(saasbb  >+IsabI‘ 


* 


m 


34 


ImtSAASBB  * 


m41  =  “m14 


m, „  =  ~m„, 
42  24 


m43  =  ”m34 


m, ,  = 


“44 


(3.63) 


In  this  case,  we  have  a  maximum  of  seven  independent 
elements.  We  note  from  (3.62)  and  (3,63)  that,  all  elements 
of  the  Mueller  matrix  [M]  are  real. 


b)  Derivation  of  the  modified  Mueller  matrix  I  Mm]  from  the 
scattering  matrix  [S]  : 

Using  the  definition  of  [Mm]  in  Eq.(3.61)  and  of  gm  in 

Eq.(3.11),  the  elements  of  the  modified  Mueller  matrix  [Mm] 
can  be  obtained  directly  from  the  scattering  matrix  [ S ( A , B ) ] 
(Appendix  F) ;  thus  for  the  bistatic  case,  we  obtain 


PAGE  59 


'saa'* 

'Vs 


IsabI‘ 

SBB  I  * 


**tSAAS*B  1 
EetS5ASBB  1 


l«SU?U  > 


2Ra(SAASBA*>  2Re<SABSBB*>  Ra<SAASBB*+SASSBA*) 

-2«saasba*>  -2I"<sabsbb*>  -i“(saasbb‘+sabsba*>  R8{smsbb*-sabsba*> 


(3.64) 


and  for  the  monostatic  case  we  only  have  seven  independent 
elements,  since 

M21  *  M12’  M31  =  2M13’  M41  =  “2M14’  M32  =  2M23’ 

M42  =  “2M24’  M43  “  "M34  and  M44  =  M33-2M12 


3.6.2  Derivation  of  the  Scattering  Matrix  [SI  from  the 

Mueller  [Ml  or  Modified  Mueller  fMml  Matrices  : 

Since  it  is  useful  to  characterize  clutter  behaviour  by 

its  optimal  polarization  properties,  it  is  desirable  to 

express  the  amplitudes  and  the  phases  of  the  scattering 

matrix  elements  S.A,  S4n,  S_.  and  S_„  in  terms  of 
AA  AB  BA  BB 

the  elements  of  the  Mueller  matrix  elements  m^  or  of  the 
modified  one  denoted  here  by  Mij . 

a )  Derivation  of  the  scattering  matrix.  [SI  from  the  Mueller 
matrix  [Ml  : 

As  shown  in  Appendix  (G) ,  for  the  bistatic  case;  the 
amplitudes  of  the  elements  of  the  matrix  [S]  are  given  in 
terms  of  the  elements  of  [M]  by  : 

|SAa'  =V1("ll+”l2'h"21+'“22) 


SAB  1 

SBA!  *\/ 4Cmll‘H'l2'm2l'"22> 


(3.65) 
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and  their  phases  by  : 

m14+m24 

*AA  =  *AB+arctan  - 

^13^23 

is  arbitrary  (“0) 


m14"ra24 


m4l‘m42 


^BA  “ 


^AB+arCtan 


+arctan 


m 


13  m23 


m3l‘“32 


*BB  “  ^AB4arCtan 


m4l"m42 


m3l"*m32 


(3.66) 


which  simplify  in  the  monostatic  case  so  that 
'V  "  |SBa'  a"d  *AB  =  *BA‘ 

b)  Derivation  of  the  Scattering  Matrix  TS1  from  the  Modified 
Mueller  Matrix  [Mm|  : 

Similary  for  the  case  of  the  modified  Mueller  matrix  [Mm] 
as  shown  in  (Appendix  H) ,  we  obtain  for  the  bistatic  case  : 


M 


14 


AA 


•‘fihl  ■  *  AA  =  ♦AE+ar>:t“ 


13 


AB  1 


YMi2  »  ^ab  is  arbitrary 


"BA1 


VM21  *  ^BA  =  *AB+arctan 


M34+M43 


M33-M44 


(3.67) 


jf.L 


•'ll  i  I  tkaift  »  >iirT  r‘  tw-i 


■  rnktSHm 
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M 


'42 


BB1 


~\'22  ’  ^BB  ""  ^An'*'arc^an 


AB 


M, 


32 


which  for  the  monostatic  case  simplifies  so  that 


SAb'  =  'V  “d  *AB  “  *BA- 


3.6.3  The  relation  between  the  Mueller  matrix  f Ml  and  the 
modified  Mueller  matrix  fMml  : 

For  completeness,  the  relationship  between  the  Mueller 
[M]  and  modified  Mueller  [Mm]  matrices  is  derived  in  this 
section.  Rewriting  (3.59)  and  (3.61)  which  defines  [M]  and 
[Mm]  matrices  as  follows  : 


SS(A,B)  =  [M]gi(A, B) 


(3.68) 


Sm£C A,B)  =  [Mm]gmi(A,B)  (3.69) 

where  g  =  (I,Q,U,V)  is  the  Stokes  vector  which  is  defined  in 
Eq.(3.9)  and  =  (i(i+Q) >i(I’Q) »U,V  )  is  the  modified 
Stokes  vector.  From  the  definitions  of  g  and 


g  ,  one  can  get  the  relation  : 
— m 


>B)  = 


i(l+Q) 

i(I-Q) 

U 

V 


i  i 

i  -i 

o  o 

0  0 


0  0 

0  0 

1  0 
0  1 


I 

Q 

u 

V 


(3.70) 


or  g  (A,B)  -  [R]g(A,B) 

-m  — 

where 


[R]  = 


i  i  o  o 
i  -i  o  o 
0  0  10 
0  0  0  1 


From  (3.71),  we  can  write  : 
SfflS(A,B)  =  [R]gS(A,B) 


(3.71) 


(3.72) 


(3.73) 
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From  (3.73)  and  (3.68),  then  : 


J  mS(A,B)  =  [R]  gS(A,B)  =  [R][M]  g^A.B) 


(R][H][R'1]gffli(A,B) 


(3.74) 


where  [R  *]  is  the  inverse  of  the  real  matrix  [R]  of 
(3.72)  which  can  be  written  as  : 


[R-1]  = 


110  0 


1-10  0 
0  0  10 

0  0  0  1 


(3.75) 


From  (3.69)  and  (3.74),  one  can  write 


-1, 


[Mm]  =  [R][M][R  ‘] 


(3.76) 


or 


[M]  =  [R  A]  [Mm]  [R] 


(3. 77- 


In  summary,  the  relation  between  [M]  and  [Mm]  is  derived 
and  is  given  by  the  pair  of  (3.76)  and  (3.77),  where  [R]  and 

[R  ]  are  given  by  (3.72)  and  (3.75),  respectively. 

3.7  Reconstruction  of  the  scattering  matrix  [SI  and  the 
Mueller  matrix  [M]  from  the  optimal  polarizations  known  on 
the  Poincare  sphere  : 

In  this  section  the  reconstruction  of  the  scattering 
matrix  [S]  and  the  Mueller  matrices  [M]  and  [Mm] is  derived 
assuming  that  the  C0P0L  null  pair  or  one  COPOL  and  one  XPOL 
polarization  null  are  given,  it  will  be  shown  that  knowing 
only  the  XPOL  null  pair  is  not  sufficient  to  reconstruct 
these  matrices  as  expected.  To  solve  this  problem, we  have  to 
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use  the  equations  in  the  previous  sections  of  this  Chapter. 
The  problem  now  is: given  two  polarization  nulls, either  two 
COPOL  nulls  or  one  COPOL  and  one  XPOL  null  but  not  two  XPOL 
nulls,  such  that  their  representation  on  the  Poincare  sphere 
are  given  by  (p , Q L ’  )  and  (p,82,02'),  where 

p=span{[S]}  is  the  radius  of  the  Poincare  sphere,  8  is  the 
colatitude  and  0'  is  the  longitude.  Then,  it  is  required  to 
reconstruct  the  scattering  matrix  [S]  with  relative  phase 
and  similarly  the  Mueller  matrices  [M]  and  [Mm] .  First,  we 
have  to  calculate  the  auxiliary  parameter  p  which  is 
mentioned  in  Section  C 3 . 5 )  from  knowing  8,  0'  at  any 
polarization  null.  Using  (3.56),  one  can  get. 


1-u 

p  =  -j  - 

1+u 

where  u  is  a  complex  number  and  it  can  be  calculated  using 
(3.57)  and  (3.58)  as  follows  : 


'  l+cos8 


u  = 


"  J  0 ' 


l-cos8 

where  (8,0*)  is  the  coordinate  of  one  of  the  polarization 
null  on  the  Poincare  sphere.  This  means  if  (8,0')  is  known 
for  such  a  polarization  null,  its  corresponding  complex 
parameter  p  can  be  calculated  using  (3.78)  and  (3.79).  Next, 
we  will  calculate  the  scattering  and  the  Mueller  matrices 
from  the  optimal  polarizations. 


3.7.1  the  COPOL  null  pair  is  known  • 

Let  the  COPOL  polarization  null  pair 

f  o  CO  ,  CON  j  /  a  CO  .  CO»  i  , 

(p»01  )  and  (p,82  ,02  )  be  given. 


(3.78) 


(3.79) 


where  we  know 


from  Section  3*5  that  p=span{[S]}.  By  using 
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(3.78)  and  (3.79)  we  can  calculate  the  corresponding  complex 

CO  CO 

parameters  p^  ,  p^  Each  of  these  parameters 

should  be  a  root  of  the  first  equation  of  (3.50).  This 
implies  for  the  monostatic  case  : 

•m+(0,#b»«»i“  sab  ■  0 

and 


Wie">,8m**<’2CO  SAB  ■  0 
where  SAA>  S^g  and  Sgg  are  the  elements  of  the 

scattering  matrix  [S(A,B)]  in  the  orthogonal  basis  (A,B). 
From  (3.80)  and  (3.81),  one  can  calculate  two  elements  of 
the  scattering  matrix  in  terus  of  the  third  one.  For 
example,  solving  both  equations  to  calculate  S^,  Sgg  in 


torms  of  S,-, 
AB 


one  can  write  : 


*AA 


0  CO  CO 
"2P1  p2 


co ,  co 
Pl  +p2 


AB 


and 


-2 


SBB 


co ,  CO 

Pl  +p2 


AB 


then  the  scattering  matrix  can  be  written  in  the  form  : 


„  co  co 
2pl  "2 


[S(A,B) ]  =  S 


AB 


co .  co 
Pl  +p2 


-2 


co ,  co 
P1  +p2 


(3.80) 


(3.81) 


(3.82) 


(3.83) 


(3.84) 
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Now  we  have  to  calculate  S^,  where  by  using  (3.52) 

|SAAi2+2,SAB|a+!SBB(a  =  P  =  span{[S(A,B)]>  (3.85) 

We  substitute  (3.82)  and  (3.83)  into  (3.85),  then 


The  absolute  phase  of  cannot  be  reconstructed  from  the 

knowledge  of  the  optimal  polarization  pairs  as  was  clearly 
\ 

shown  by  Kennaugh(1952) .  The  scattering  matrix  with  relative 
phase  only  can  be  reconstructed. Letting  then 

S^g=|S^g|.  Substituting  from  (3.86)  into  (3.84),  one  can 

write  the  reconstructed  scattering  matrix  [S]  with  relative 

phase  in  terms  of  the  COPOL  null  pair  (p^co,p2co)  as 

follows  : 

[S (A, B) ]  =  K 

where 

K=fi/1{ ( P  XC°+P2C0 1 a+2 1 p 1cop 2C° | *42 } ' , 

a  =  -2p1c°p2coexp(-j0E), 

b  =  |pj  +P2  I , 
c  =  -2exp(- j0g)  and 

$E  -  phase(p1  +p2  ) 


— 

a  b 
b  c 


3.7.2  One  COPOL  and  one  XPQL  null  are  known  : 


CO  X. 


Let  (p^  ,p^  )  V  xnown  on  the  Poincare  sphere  of 
radius  p,  where  p^r°  is  aay  COPOL  null,  i=l  or  2  and 

X  CO 

pj  is  any  XPOL  null,  j~l  or  2.  p^  should  satisfy 


the  first  equation  of  (3.50)  and  p^  satisfy  the  second 
one  T1  '  .  .eans  : 


and 


SAA+('iC°),SBB+2'’lCO  SAB  *  0 


mulitply  (3.33)  by  (Pj  )  and  a^d  it  to  (3.89),  then  : 

„  co>  a..*  |  Xij,. 

2pi  (p^  )  -IPj  r+i 


BB 


AB 


(pp'Cp/'V+Pj* 

Substitute  (3.90)  into  (3.88),  then  : 


CO  CO  |  Xi]  ,  x 
- a  .  I o .  !  -2p 


Pi  - p±  iPj 


.1 


SAA  “  SAB 


co 


,  X  *,  CO. 2,  J 
(Pj  )  Cp4  )  +Pj 


Using  (3.90)  and  (3.31)  nnd  (3.85),  then  we  can  write  : 

|SABi  =  P/D  <  |  (pjX/(piCO)a-»PjX|  > 

wi.ere 

D  *  {2|(pjX)%iC°)2+PJX|* 

+  lpiCO|*llrico-olCO|pJX|«-2pj*H* 

+  l[2piC°(p.V-|pjX|2+l]|a) 

From  (3.90),  (3. 91)  and  (3.92),  we  can  write  the  scattering 
matrix  with  relative  phase  C^^g=D)  as  following  : 


(3.88) 


(3.89) 


(3.90) 


(3.91) 


(3.92) 


(3.93) 
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[S(A,B) ]  =/  - 

D 


where 


a  b 


c  d 
.  j 


D  is  given  by  (3.92), 
a  =  PiC°[PiCO-piCOipjX|i-2pjx]e-^E 

b=  I (pjX)*(p1C°)a+PjX| 

o=  I (PjX)*(PiC°)l+PjX| 

d  =  '[2n.co(p,x)*-|Pjx|*+i]e"^E 

i 

*E  =  phase  of  [ (PjX)*(piC°)J+pjX] 

3.7.3  The  XPOL  null  pair  is  known  : 

X  X 

Let  (p^  ,p£  )  be  the  XPOL  null  pair,  then  each 
value  should  satisfy  the  second  equation  in  (3.50),  where 

and 

-("2X)*V'2XSBBJ-SABt1-'!2X^2X)*)  -  0 

Multiply  (3.97)  by  p^  and  subtract  the  result  from 
(3.96)  after  multiplication  by  p2x,  then  : 

«AbI^18^2*^2X5%iX+P2*-Pi^2X<«»1*5*]  -  0 

Using  the  relation  Pj*^*)*  =CPiX)*P2X=  “1 

i.e.  the  property  of  the  XPOL  nulls  to  be  antipodal, 


(3.94) 


(6.95) 


(3.96) 


(3.97) 


(3.98) 


i*ie  i  m"  fcU 


-  b- in.i.  ■> 


mmWhhhm.  hi 


which  can  be  derived  easily  from  (3.55a)  and  (3.55c). 
£q.(3.98)  will  vanish  for  any  value  of  p^x  and  p^X. 

This  is  also  true  when  one  tries  to  calculate  S__  in  term 

DD 

of  This  means  the  scattering  matrix  [S]  cannot  be 

determined  by  using  the  two  XPOL  nulls  as  was  previously 
established  by  Kennaugh  and  shown  also  in  Huynen  (1970). 

In  summary,  the  scattering  matrix  [S]  with  relative  phase 
at  fixed  aspect  and  at  a  given  frequency  can  be 
reconstructed  by  using  (3.87)  for  the  case  in  which  the  two 
COPOL  nulls  are  known,  and  by  using  (3.94)  for  the  case  in 
which  one  COPOL  and  one  XPOL  null  are  known.  We  note  here 
that  [S]  cannot  be  reconstructed  if  only  the  two  XPOL  nulls 
are  known.  After  calculating  the  scattering  matrix  [S],  one 
can  calculate  the  Mueller  [M]  and  the  modified  Mueller  [Mm] 
matrices  in  the  monostatic  case  using  (3.63)  and  (3.64), 
respectively. 
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CHAPTER  FOUR 
NUMERICAL  RESULTS 

4.1  Introduction  : 

In  this  Chapter,  we  will  use  the  theory  of  Chapter  Three 
to  calculate  the  two  unknown  matrices  given  either  the 
scattering  [S],  the  Mueller  [M] ,  or  the  modified  Mueller  [Mm] 
matrices  for  different  targets  and  for  sea  clutter.  We  will 
concentrate  on  the  monostatic  relative  phase  case  for  the  [S] 
matrix.  It  should  be  noted  that,  the  bistatic  case  also  can 
be  calculated  by  using  the  formulas  of  the  same  Chapter. 
Also,  the  CQPOL  and  XPOL  nulls,  at  a  given  aspect  and  fixed 
frequency,  and  their  representation  on  the  Poincare  sphere 
for  each  case  are  calculated.  It  should  be  noted  that,  for 
calculating  the  scattering  matrix  [S]  from  the  Mueller 
matrices  [M]  or  [Mm]  in  the  monostatic  relative  phase  case, 
only  seven  elements  from  [M]  or  [Mm]  are  needed,  e.g.  m^, 

m12*  m13*  “l4’  “22*  m23  and  m24  aS  was  shown  in 

(3.63)  and  (3.64). 


4.2  Targets  with  Simple  Shapes  :  (Huynen  1970) 

Example(l)  :  Largs  metallic  sphere  or  a  flat  plate  : 

In  this  example  a  large  ideally  conducting  sphere  or  flat 
plate  at  normal  incidence  is  considered.  The  target  shape  is 
shown  in  Fig.4.1.  The  scattering  matrix  is  given  by  : 


li;  j  1 1  <1j ul'h J.'tilt ,'lW l a! uk kllii ..  , 


'IjLA./il/t.ihU  y/kw 
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the  Mueller  matrices  are  given  by  : 

_1  0  0  O' 

[M]  =  [Mm]  *0100 
0  0  10 
.0  0  0  1. 

Optimal  polarization  :  (t’ig.4.2) 

a.  COPQL  nulls  : 

co  .  co  .  o  .  co  , 

P1  3  j,  3  0  ,  arb. 

CO  ■  a  GO  i  a  O  j  GO  « 

P2  3  -J,  02  180  *  *2  arb‘ 

This  means  the  COPOL  nulls  lie  at  the  North  and  South  poles. 

b.  XPQL  nulls  : 

They  exist  anywhere  at  antipodal  locationsi  on  the  Equator 
(9=90°  major  circle). 

In  this  example,  if  the  polarization  of  the  incident  wave 
is  loft  circular  then  the  return  signal  will  be  right 
circular  and  vice  versa.  This  means  if  both  the  radar 
transmitting  and  receiving  antennas  are  adjusted  to  use 
left(or  right)  circular  polarization,  then  t\a  receiver  will 
receive  no  (or  minimum)  return  signal  from  the  target.  Also, 
the  received  return  signal  will  be  maximum  if  both  antennas 
are  uaing  the  same  linear  polarization. 
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Example  2:  Metallic  trough:  (Fig. 4. 3) 

The  target  in  this  case  is  &  large  metallic  trough  (two 
planes  intersecting  at  90°)  oriented  with  axis  (the  plane's 
line  of  intersection)  horizontal  or  vertical.  The  view  angle 
is  considered  normal  to  the  trough's  open  surface.  This 
target  has  a  two-bounce  reflection  characteristic.  The 


scattering  matrix  is  given  by  (Huynen  1970)  : 


and  the  Mueller  matrices  are  : 


[M]  =  [Mm]  = 


10  0  0 
0  10  0 
0  0-10 
0  0  0  -1 


Optimal  polarizations:  (Fig.  4.4) 
a.  C0P0L  nulls  : 


=  -1, 

0^°  = 

90°, 

<t>xCO  •  -90 

=  1, 

A  CO 

02  - 

90°, 

=  90° 

b.  XPOL  nulls  : 

They  exist  anywhere  at  antipodal  locations  on  the  major 
circle  $  =  0°. 

NOTE  :  In  comparision  with  the  sphere  results  of  Example  1, 

we  note  the  important  difference,  in  phase  of  element  S,.„ 

*3  0 

for  Example  2;  i.e.;  the  relative  phase  relating  the  two  co- 
polarizsd  elements  is  of  paramount  importance. 
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Example  5  :  Linear  target  oriented  by  an  angle  j*  to  the 

horizontal  (x-axls)  (Huynen  1970)  : 


The  target  model  is  shown  in  Fig. 4. 9.  The  scattering 

matrix  is  given  in  terms  of  the  angle  $  by  : 

cos2V  sin^cosi/*] 
sini/»cos^  sin2^ 

The  Mueller  matrices  are  calculated  using  (3.63)  and  (3.64) 

by  : 


[S«0] 


[M«»)l  «  i 


1  cos  2  ip  sin2t}» 

cos2t{>  cos22i/i  £sin4</> 
sin2  4>  isinAip  sin22if» 


[Mm(i/»)]  " 


COS  *1^ 

£sin22ip 

2cosIi£sin^ 

0 


isin22tj> 

sin4^» 

2sin>tpcost|) 

0 


cos,<psin4» 

sinJV»cos# 

isin22if> 

0 


Optimal  polarizations  : 

a)  C0P0L  nulls  : 

PlCO=P2CO=-cot*,  eiCO-e2CO-90#,  01CO=02CO=2^±1T 

b)  XPOL  nulls  : 

Plx=tan*,  p2x=-cot^,  01x=02x=9O°>  **=*^=2+ 


Special  cases  : 

(1)  The  target  is  aligned  horizontally  (i/>g0)  ; 


[S(*-0)J 


1  0 
0  0 


PlCO=P2CO=--,  01CO=82CO=9O°>  ^c%2co=ir 

PlX-0.  P2XS— ,  01X302Xa9Ot’,  ^X^2X»0 


The  Co-Pol  and  X-Pol  nulls  are  shown  in  Fig. 4. 10. 
(ii)  The  target  is  aligned  vertically  (^=90°)  ; 


(S(*»90°)] 


r  « 
0  0 

0  1 


Fig.  4.12  :  COPQL  and  XPOL  nulls  for  a  linear  target  with'f'=45  . 
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4.3  Targets  for  more  Complex  Shapes  :  (Crispin  et  al  1961) 
Three  targets  for  more  complex  shapes  are  taken  into 
consideration  in  this  section.  The  data  are  taken  from 
(Crispin  et  al  1961)  for  three  target  models  :  missile  with 
and  without  fins,  and  the  ncse  cone  model.  The  measured  data 
for  the  radar  cross-section  and  relative  phases  are  shor-m  in 
Table(4-1).  The  radar  cross-section  were  measured  relative  to 
a  conducting  sphere.  The  geometry  for  the  cross-section  study 
used  in  The  University  of  Michigan,  Radiation  Lab.  Report 
(Crispin  1961)  is  shown  in  Fig. 4. 13. 


Target 

Aspect 

RCS  in  DBS  Relative  to 

Sphere _ _ _ _ 

Phases  in 
Re  1 ,  to  A 

Sphere 

Diameter 

a 

XX 

0 

xy 

a 

yy 

X 

X 

-e- 

'  yy 

Nose 

Cone 

26° 

3.6 

-16.1 

2.8 

128,0° 

121.0° 

2.945" 

Missi  le 
with 

FI  ns  ' 

o 

O 

m 

-9.3 

-20.5 

i 

O 

CO 

98.60° 

159.3° 

II 

1 .93 

Missile 

without 

Fins 

o 

O 

CO 

5.9  ‘ 

-30.0 

6.3 

34.5° 

35.0° 

II 

1 .98 

Table  4 • 1 :  Experimental  data  for  RCS  (Crispin  et  al,  1 9^ 1 ) 


The  frequency  is  9.7  GHz  and  the  range  R=36-£ft.  for  the 
nose  cone  target  model  and  R=33ft.  for  the  missile  model. 
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Fig.  4.13  :  Penmetry  for  cross  section  study  (Crispin  et  al  1961). 
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Example  (8):  Missile  model  without  fins  ; 

The  geometry  of  this  model  is  the  same  as  the  last  one  but 
without  fins  (Fig. 4. 16).  The  scattering  matrix  is  calculated 
from  the  data  given  in  Table  4.1  at  aspect  angle  180°.  The 


calculated  scattering  matrix  relative  to  a  metallic  sphere 
with  diameter  1.98"  is  given  by  : 

-1 


0.111x10 


0 . 569+j0 . 391 

0 . lllxlO*1  0 . 592+jO . 415 


[S]  - 

and  the  Mueller  matrices  are  : 


[M]  = 


0.5 

-0.23x10 


-1 


0.129x10 

0.261x10' 


-1 


-0.23x10 

0.5 

-0.257x10' 
-0 . 892xlo" 


0.  129x10 
-0.257x10 
0.499 
0.436x10 


-1 

-3 


-2 


-0.261x10 
0.892x10 
-0.436x10 
0.499 


-3 

-2 

-2 


[Mm]  * 


0.477 

0.123x10 

0.126x10 


-3 

-1 

,-2 


0.123x10 

0.523 

0.131x10 


-3 


-2 


[-0.866x10  0.918x10 

Optimal  polarizations  :  (Fig. 4. 18) 
a)  COPOL  nulls  : 


-1 

-2 


0.63x10 

0.656x10' 

0.499 

0.436x10' 


0.433x10 

-0.459x10 

-0.436x10 

0.499 


-2 

-2 

-2 


co 


■1 


p. ''“=-0.1669x10  *- jO  . 9684,  8 . CO=177 . 921°,  '*‘‘'*-28 . 335 


co 


P2CO=-0. 839x10  N-.i0.986,  B^^O .  934°  ,  *  ww*-31 . 12' 

b)  XPOL  nulls  : 

0:X=9O.568°,  ^*=150.796° 


.-1. 


co. 


CO. 


Px  =3. 837-j0. 078, 


P2X=-0.2605+j0.5298,  92X=89.431°,  $2X=-29.204° 


'ill  t  ■  'imri  ** 
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4.4  Sea  clutter  : 

Example  (9) :  Sea  clutter  and  noise  :  (Fig. 4. 19) 

The  data  of  this  example  is  taken  from  (Daley  1978)  in 
which  he  based  his  calculations  for  the  sea  clutter 
scattering  matrix  on  (Valenzuela  1968) .  The  parameters  in 


this  case  are  given  by  : 

a)  dielectric  constant  of  the  sea  water  =70-j70 

b)  wind  velocity  =10  m/sec 

c)  tilt  angle  (for  slightly  rough  surface)  =4° 

d)  depression  angle  =10° 

e)  aspect  angle  =120° 

f)  propagation  constant  =68.30° 

g)  noise  level  =-89.1db 

the  scattering  matrix  is  given  by  : 


[S] 


0.347xl0“1+j0.206xl0'1 
0 . 686x10 " X+j 0 . 737xl0~3 


0. 686x10“  VjO. 737x10 
0 . 951+j0 . 291 


where  according  to  Section  3.6.1 


[M]  = 


and 


[Mm]  = 


Optimal  polarizations 
a)  COPOL  nulls  ; 


0.5 

-0.494 

0.679xl0_1 

-0.179X10'1 

■0.494 

0.491 

-0.631X10"1 

0.206X10'1 

0 . 679xl0_1 

-0.631X10'1 

0 . 436xl0~ 1 

0.946xl0-2 

0.179X10"1 

-0 . 206x10” 1 

-0 . 946x.lo”2 

0. 342x10“ 1 

0.162xl0“2 

0.47lxl0“2 

0 . 239xlQ-2 

0 . 138xl0~2 

0.471xl0"2 

0.989 

0.655X10'1 

-0.193X10'1 

0.479xl0"2 

0.131 

0.436xl0_1 

0.946xl0“2 

-0 . 277xl0“2 

0 . 385xl0_1 

-0.946xl0'2 

0. 342x10 “1 

(Fig. 4. 20) 


p . C°=-0 . 3436* 10*1- jO . 1713 , 


0,C°=lO9.414c 


<hCO=-b.0S5* 


p„CO=-0 . 9805x lO*1- jO . 2102 , 


02CC=66.486°, 


*2CO=-11.708o 
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b)  XPOL  nulla  : 

P1x»i3.6727-j3.6on,  8  x«92.054°,  0^172.172 

p^-O.esapxlO^+jO.iaOlxlo'1,  82X»37.946#,  12X*-7.828# 

Example  (10):  Sea  clutter,  target  and  no is a  ; 

In  this  example,  the  conditions  of  the  sea  clutter  and 
system  noise  are  the  same  as  in  the  example(9)  but  a 
simulated  target  is  present  and  no  target-clutter  interaction 
is  assumed.  The  scattering  cross-sections  of  the  simulated 
target  arc  given  by  (Daley  1978)  : 
a^h(target)  *  o^(clutter) 

a  (target)  *  a, ,  -3db 
w  oh 

ahv(target)  *  <3hh‘10db 

and  the  phases  are  calculated  randomly. 


The  total  scattering  matrix  in  this  case  for  the  same 
parameters  of  example  (9)  is  given  by  : 


-1  -3l 

0. 161+j0.69xin  0. 257- jO. 171x10 

[S]  = 

0. 257-jO. 17:xl0"3  0 . 817+j0 . 412 

The  Mueller 

matrix  is  : 

"o.5  -0.403  0.251 

-0. 884x10 -1 

[Ml  - 

-0.403  0.368  -0,169 

0.124 

0.251  -0.169  0.226 

-0 . 99xl0*2 

-0. 884x10* 1  -0.124  0.99xl0*2 

0.935xlQ~1m 

The  modified  Mueller  matrix  [Mm]  is  given  by  : 

r 0.306x10  1 

0.662X10*1 

0.413xl0-1 

0.178x10* 

>  -1 
662x10 

0.837 

0.210 

-0.106 

1  . . 826x10 

0.420 

0.226 

-0.990x10* 

1  m  1 

1.-0.355x10 

0.212 

0.990xl0*2 

0.935x10* 
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Optimal  polarizations  (Fig.  4.21)  : 
a)  COPOL  nulla  : 

PlCOa-0. 1836-jO. 2583,  01co«118 . 001° ,  4*°=- 21.209° 

P2CO—0.3182+j0.S121,  92CO«41.3Ua,  ♦2°°«-4S  .001° 


b)  XPOL  nulls  ; 

PlX*3. 1385- j  1.1042,  6  ^“lOQ.  543°,  *^148.062* 

P2X=-0 . 2835+jO .9975x10*^,  82*=79 .457* ,  *2X»-31.938° 


la  the  above  two  examples , there  is  a  difference  in  two 
elements  of  the  modified  Mueller  matrix,  e.g.  and. 

as  compared  to  Daley's  results  due  to  a  printing  mistake  in 
Ishimaru's  book  (1978,  vol.l,  pp.35)  and  the  correct  elements 

should  according  to  Appendix  (F)  ^3  i*2Re^  saaSBA  ^  ’ 

M^^-ZImiS^S^  }.  In  Daley's  results,  he  used 
★ 

instead  of  S.,  .  Also,  it  should  be  noted  that  some  of 
BA 

the  elements  of  [Mm]  differ  in  sign  with  the  ones  mentioned 
in  (Ishimaru  1978)  and  (Daley  1978)  due  to  the  difference  in 
sign  of  g^  la  Eq.(3.9).  Also,  we  noticed  that,  the  XPOL 

nulls  are  antipodal  and  they  should  lie  on  a  major  circle 
with  the  COPOL  nulla  and  bisect  the  arc  between  them.  But  in 
our  model  example  the  XPOL  nulls  are  shifed  with  very  small 
angles  in  both  example  which  may  be  due  to  the  fact  that  the 
scattering  matrix  does  not  completely  satisfy  the  relative 
phase  concept  or  due  to  some  measurement  errors  which 
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requires  further  analysis.  For  example,  the  scattering  matrix 
in  the  first  example  has  0.62°  absolute  phase  and  for  the 
second  one  0.04°.  We  emphasize  here  the  potential  use  of 
calculating  the  C0P0L  and  XPOL  nulls  from  measurement  data 
immediately  during  measurement  campaigns  for  the  purpose  of 
checking  the  accuracy  of  the  measurements.  This  aspect  is 
further  discussed  in  Section  4.5;  and  will  be  treated  in  all 
detail  in  a  forthcoming  report. 


Fig.  4.21  :  C0P0L  and  XPOL  nulls  for  a  simulated  target  and  sea 
clutter. 
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4.5  Discussions  ; 

The  numerical  investigation  of  some  target  shapes  and  sea 
clutter  has  shown  the  feasibility  of  applying  the  theory  of 
Chapter  Three  to  calculate  the  two  unknown  matrices  of  the 
three  matrices  [S],  [M]  and  [Mm]  if  one  of  them  is  known.  We 
concentrated  on  the  aonostatic  case  only  which  is  of  specific 
interest  here.  It  should  be  noted  that,  the  scattering  matrix 
with  relative  '  phase  can  be  reconstructed  from  [M]  or  [Mm] . 
There  are  seven  elements  required  for  calculating  [S]  which 
«.  mu,  b12>  .lv  aiv  m22.  «23  and  «24- 


For  calculating  the  optimal  polarisation  in  the  previous 
examples  it  is  noticed  that,  as  expected,  the  cross* 
polarization  (XPOL)  nulls  are  antipodal  on  the  Poincare 
sphere  and  they  bisect  the  angles  between  the  co-polarization 
(CQPOL)  nulls.  Also,  it  should  be  noted  that  [M]  and  [Mm] 
have  only  seven  independent  elements  as  shown  in  the 
calculations  for  the  monostatic  relative  phase  case  as  shown 
in  Section  3.6.1. 


More  studies  are  needed  for  analyzing  the  same  shapes  of 
targets  and  others,  and  also  for  sea  clutter  with  and  without 
target  to  see  the  effect  of  changing  the  frequency  or  aspect 
angle  on  the  COPOL  and  ^?0L  null  locations  and  how  they  are 
moving  on  the  Poincare  sphere.  For  the  mere  reason  to  keep 
thin  report  of  still  modest  size,  aspects  of  target  optimal 
polarization  characteristics  will  be  treated  in  detail  in  a 
forthcoming  report. 
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CHAPTER  FIVE 

CONCLUSIONS  AND  RECOMMENDATIONS 
V*  have  thoroughly  reviewed  the  literature  on  the  basic 
theory  of  polarisation  utilisation  in  radar  target 
allocation,  dataction,  imaging  and  idantification.  Basad  on 
this  study,  we  dacidad  to  davalop  tha  thaoriaa  and  computer 
ass is tad  numarical  algorithms,  stap  by  stap,  rachacking 
every  altaraata  raprasantation  found  in  tha  litaratura  and 
thus  correcting  several  misrepresented  f emulations .  We  will 
retain  this  approach  rhr*.  nur  studies  and  produce  four 

major  interim  reports  annually. 

3.1  Progress  Reported  in  this  Report (January  1j,  1981) 

Wa  rederived  tha  basic  formulations  of  polarisation 
raprasantation  for  tha  radar  case.  Wa  established  tha 
relationship  between  tha  2x2  radar  scattering  matrix[S],  tha 
4x4  Stokes  reflection  matrix[M],  and  tha  4x4  modified 
Mueller  matrix[Mm] ,  and  vice  versa  for  both  tha  monostatiu 
and  tha  bistatic  casas.  Wa  derived  tha  associated 
expressions  for  calculating  tha  CO -POLARIZATION  (COPOL)  as 
wall  as  CROSS -POLARIZATION  (XPOL)  nulls  and  their 
presentations  on  tha  Poincare  sphere.  Inversely,  given  tha 
radius  p*span{[S]}  of  tha  Poincare  sphere  and  tha 
coordinates  of  either  both  COPOL  nulls,  or,  of  one  COPOL  and 
one  XPOL  null,  we  derived  the  expressions  for  tha  elements 
of  [ S ] ,  [M]  and/or  [Mm]  re-expressed  as  functions  of  p  and 

the  coordinates  of  the  two  respectively  given  nulls. 


PACffi  98 


We  have  developed  computer  program  algorithma  for 
calculating  thasa  scattering  matrices,  the  associated 
optimal  polarisations  '.’or  the  monostatic  and  bistatic  cases 
for  a  variety  of  perfectly  conducting  target  shapes.  We  have 
initiated  close  collaboration  with  several  research 
laboratories  in  possession  of  excellent  measurement 
facilities  for  obtaining  reliable  monostatic  as  well  as 
bistatic  scattering  data,  and  some  of  the  data  made 
available  have  been  used  to  calculate  the  respective  optimal 
polarizations. 

We  have  extended  the  mono-chromatic  theory  of  the  optimal 
polarization  concept  to  the  quasi -coherent  and  pen-chromatic 
cases,  and  ws  are  completing  a  study  on  a  novel  vector 
scattering  approach  of  extracting  the  useful  target  signal 
from  cluttar  perturbed  data  which  will  be  presented  in 
detail  in  one  of  the  forthcoming  reports. 


5.2  Conclusions 

In  analyzing  the  model  data  used  for  verification  of  She 
optimal  polarization  concept,  it  becomes  very  evident,  that 
we  require  to  measure  in  the  monostatic  case  the  relative 
phases  and  magnitudes  of  the  elements  of  the  scattering 
matrix  [S]  to  describe  the  complete  electromagnetic 
properties  of  a  target  uniquely.  For  example  (see  Section 
4.2),  if  the  relative  phase  between  the  two  co-polarizsd 
components,  and  Sgg,  is  not  known,  we  find  that  for 


grossly  different  shapes  the  magnitudes  can  be  identical 
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resulting  la  ambiguity  of  identification  (Section  4.2, 
examples  1,  2  and  S).  Thus,  the  relative  phases  existing 
between  all  elements  of  the  scattering  matrix  t s 1 shK  ue*^ 

to  be  measured  for  proper  unique  presentation  of  the  target 
polarisation  characteristics  which  is  being  clearly  verified 
by  the  examples  presented  in  Section  4.2. 

Ve  also  note  that  the  inverse  process  of  calculating  the 
elements  of  [S],  [M]  and/or  [Mm]  from  the  coordinates  of  the 
respective  optimal  polarization  nulls,  provides  deep  inside 
into  target  characteristic  properties  and  we  should  have  a 
powerful  tool  for  introducing  the  concept  of  target 
polarization  synthesis.  In  target  polarization  synthesis, 
we  would  wish  to  design  the  shape  and  properties  of  a  target 
such  that  it  produces  a  given  set  of  optimal  polarization 
pairs. 

5.3  Recommendations 

Based  upon  the  results  of  our  model  verification  studies, 
we  conclude  that  it  is  meritorious  to  further  advance  the 
optimal  polarization  concept  first  introduced  by  Professor 
Edward  M.  Kennaugh,  subsequently  extended  by  Dr.  J.  Huynen, 
and  to  utilize  its  great  potential  in  radar  target  imaging 
and  idantif  ication. 

We  also  recommend  that  in  microwave  remote  sensing  using 
both  passive  (radiometry)  and  active  (SAR,  SLAR, 
SCATTEROMETRY')  methods,  measurements  of  the  complete 
monostatic.,  relative  phase  scatter ing  matrix[S]  are  made, 
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i.a.  both  wnplitudaa  and  ,  the  relative  phases  of  all  four 
(thru*')  elements  need  to  be  executed  to  obtain  unique 
information  on  tha  scattering  properties  of  remotely  sensed 
targets.  Tha  associated  question*  of  applicable  measurement 
techniques  will  ba  tract  ad  in  a  forthcoming  :  -sport  in  grant 


detail . 
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Appendix  (A) 
Proof  of  Eqs . (3-8) 
From  Eqs. (3.3)  and  (3.7),  wa  hav«  : 
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Appendix  (B1) 

Stokaa  Paraaatars 
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«2  *  2Ra(hxhy*>  ■  0  (B-l) 

g3  •  -2Ina{hxhy*}  «  V 

From  Eqa.(A-4),  (A-S),  (A-6)  and  (A-8)  in  Appandix  (A),  wa  can 
rawrlta  tha  Stokes  paranatars  in  tarn*  oil  (a  ,a  ,5  ,6  ) 

3a  X  y 

or  (a,A,t,aO  as  following  : 

«q  -  |hxl*  +  Ihyl1  •  axa  +  Ay*  -  a1  -  I 

*1  “  IV*  "  lhyl*  "  ***  "  4y*  “  *acoa2tcoa20  *  Q 

S2  «  2Ha{hxhy*}  ■  2axayco«8  »  a*coa2Tsin2*  *  0  (B-2) 

*a  *  "  24xAy*ia6  “  Al#ia2T  “  V 

Free  (B-2),  wa  can  writs  : 

gi^+Sj1^1  M  Qa  +  (Ja  *•  va 

*  44{co»12tcoa12#+coa*2Tsina2$+sinl2t} 

*  a‘‘{co»12t(coaa2ijHTiin*2A)-HiinI2t} 

■  A*(coaa2'C"Htin12T} 

*  a*  *  gQa  *  I1 

TLon 

*0*  *  *l*4‘f2*4'83t  *  *  Qa,Hla+V*  (B-3) 


■  Mm  I  ILu  'WkIJ.  ( W  1'i 
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Given 

&  -  W  +  hyhy  (C  •: 

where  h  ,  h  ere  .a  general  coop  lex  numbers,  end  h 

*  i 

end  h  ere  the  unit  vectors  in  &  end  7  directions, 

^  . 

respectively.  If  we  introduce  the  auxiliary  couple* 
parameter  u  using  h^  end  hy  components  of  (Ol)  then  : 

VJhy  l"Jp 

u  -  — . —  - -  (o; 

V^hy  1+^P 

where 


p  m  — — . .  m  polarization  ratio  (03 

hx 

But  we  have  from  Eq.(3-3)  : 

h  ■  a  *^3x,  h  *  a  e^3y  (o^ 

xx  y  y 

where  a  .a  .5  .5  are  real  numbers, 
x  y  a  / 

Then  : 


where  5  *  3  -3 

y  x 


Substitute  (05)  into  (02),  then  : 


1-JP  l-jUy/a^3 

u  ■  - - -  — — . 

1+jP  l+j(*y/*x)«J4 


CC-5 

(C-3 
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I  -  j  ( dy/ {  cos  5+ j  a  iafi  > 
l+j ( *y/®x) ( co» 5+ j » ia5 > 

Tharafora 

[l+(«y/4x)siflS] -j («y/*x)co«5 

U  "  - - - - - - 

I I  -  ( 4X)  s  lafi  ]  4- j  ( *x )  co«  5 

Introduca  tha  r*al  nuobars  X^,X^  and  Y  such  that  : 

*  l+(dy/ax)3inS 
X2  -  l-Ca^/a^siud 
Y  «  (dy/ax)co«a 
Th  an  : 

VJY 

u  ■  "  -  —  ■■" 

X24jY 

Fro*  (C-9)  : 

X^+Y*  <l+(ay/ax)sin6>*+<  (a^/a^cosaj* 

|u|  *  * - -  *  - . . . . —  ■  ■  "  . . .  ■ 

X2*+Ya  {l-Uy/a^sinfi^-K  (ay/ax)cosa}a 

1+2  ( a^  ax)  s  ina+(  ax)* 

m  ■■■■■■■  ■■'  — .  ■  ■  ■■■■  ■■■■■■■■■■■■ 

l-2(ay/ax)sina+(ay/ax)a 

Tharafor*  : 

|u| *-l  {l+2(ay/ax)3ina+(ay/ax)*>-{l-2(ay/ax)sina+(ay/ax)*> 

{l+2(ay/ax)3ina+(ay/ax),>-{l-2(4y/ax)3iaa+(ay/ax) 
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Mdy/a^sinfi 

2[l+(yax)1] 


2«  a  siafi 
x  7 


2a  a  ainfi 
x  7 


From  Eqa.(B-2)  Appendix  3  : 

g_  »  A1  ■  A  3+A  * 

*0  x  7 

g.  *  a2coa2tcoa2$  *  a  2  -  a  2 
i  *  y 

gj  ■  a2coa2tain2A  *  2x^^0035 


g,  ■  a2sin2t  ■  2a  a  sin6 
j  x  y 


Than  : 


g-  a2sin2t  2a  a  sinfi 
»  *  y 


*0 


2a  a  sin 5 
*  7 


or  ain2t  * 


From  (C-10)  and  (C-12) ,  therefore 
M2-l 

ain2t  *  — - - 

|u|2+l 

If  8  *»  it  -  2t,  therefor# , 

|u| 2-l 

coafl  *  cosC^-ar)  *  sin2t  »  ■'  -- 

M  2+i 


t 
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Re(u) 
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Proof  of  :  (I)  Dat{(S  (A  ,B  )]}  -  Det{ [S(A,B) ] >  *  invariant 

(II)  Span{ [S' (A* ,B' ) ]}  ■  Span{ [S(A,B) ] }  *  invariant 
w*  hava  from  (3.50)  ? 


*s*l 


*’*'*•  ’  <1+f>  )  tSu+P*SM+p(Su+8u)l 


s'a-b'  ■  u+p/)"1t-p*sM^saB+sAJ-pp’'sBAi 

S'B,A,  »  Cl+p(i*)'l[-P*SM+pSM+SB4-pp*SiB) 

s’b'b'  *  (i+pp*)'1tp*'sAB+sBB-P*csli+s8A)i 

Proof  of  (I)  :  D«t{ fS' (A1 .B' ) 1 )  *  D«t( fS(A.B) 1 >  *  invariant  : 


D«tC[s'(A\a')i>  -  s'A,A,s'B,B,  -  s'a,A.s'A,B, 


CI+PC  )  '■nSuVSBBMSAB+SBA)][p  CS^+S^)] 


-t-p*SM-PpSBB+SAB-Pp*SBAn-P*SAA+pSBB+SBA-p(,*SAB]} 

■  a+pp"')‘2(p*,sAA,+p,sBB,-pp*(sAa+sBA)1 

+(l.+p"ilP!)SAASBB+p*SAA(SAB+sBA)  (pp*-l)+pSM(SAB+SBA)  (l-p*p) 

•i,*!sm,+,|,1saasbb+,,*s4asba'pi’*,smsab 

+p>**SMSBB'pISBB,',!SB45BB+|,,|’*SABSaB 

+;*SiASA3'pSAB3BB'3ABSBA+|lp*SAB‘ 


‘PP  lSAASn,+P7'P  StJAST111-p*p#iS.tlS11.+pp  S  *} 
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-  (l+pp*)'2{SAASBB(l«h»1p*l+2pp*)+SABSBA(-2pp*-l-p*p**) 

★  *  ★,  *  *, 
+sabsaac“p  pp  +p  )+sbasaa('p  +PP  2+p  -pp  2) 

+SABSBB  1+p+p  2 p*‘p  5  +SBASBB  ^  "P#P  *+P  ~P+P  *  P*>  * 


(l+pp*)“2{SAAS„(l+ppw)l-SARSTlA(l+pp")1} 


AA  BB 


AB  BA 


*  (SAASBB*SABSBA* 

-  d«tt{[S(AtB)]> 

Tharafora,  dat{ [S' (A' ,B' ) ] }  *  dat( [S(A,B) ] }  *  invariant 


Proof  of  II  :  Scanas'CA'  .B')1)  »  Scant fS(A.B) 1 }  »  invariant 
Fr tM  tha  dafinition  : 

Sp«n  ( [S' (A' ,B 1 ) 1 }  *  IS'a,a,|2+|S'a,b,|x+|S,b,a,|2+|S,b,b,|1 

*  s* a'a'^'a'a*  ^'a’b^'a'b'  ^'b'a^'b'a'  ^Vb'^b'b' 

Usa  Eqs.(D-i)  to  (D-4)  and  thair  conjugatas  into  Eq.(D-5), 
tharafora  : 


★  ★.  ★  ★ .  ★  * 


Spaa{  [S'  (A*  ,B' )  ] }  -  (  IS^VS^+pCS^+S^)  ]  [S^  +o  2SBB  +p  (S^  +S]JA  )] 


+(’P#SAA'HpSBB4,SAB"pp*SBA)C’pSAA*+p*SBB*+SAB*’pp1,fsBA1'r) 
+<“P#SAA+PSBB+SBA"PP*SAB)C'PSAA*+P*SBB*+SBA*'PP  SAB  } 


★  ★ 


*  * 


+tp  ‘SAA+SBB-"  <SAB+SBA)“,,‘SAA  +SBB  '^AB  +SBA  ’X'1*'"’  > 


Tharafora, 


Cl+PP  )2span{[S,(A',B')]> 


*  ★  ★  ,  •k'M 

SAASAA  ^1+pp  +pp  +p  P  J+S™S*n  tP  P  +pP  "^PP  +1J 


BB  BB 


+'3ABSAB*^pp*+1+plp  l+PP*J+SBASBA  ^PP  +p*p  *+1+pp  1 


+SAASAB  ^p  “p  +pp  *"pp  *1+SbbS*o  tP*P  +P"P*P  *Pl 
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+SAA~SAB Ip-p+pV-pV]  +Sbb\b (pp*’V-pp**-  * 


BB  AB1 


+SA4SBA*I|,**'P*1-p‘-PP*,]+SBBSB4*[p»p*.p*p%.p; 

+Sm‘SBa'  *'*-P  -P*P 1 1  +SBB*SBAt  PP*1  -p’W-p*] 

+SAA*SBBt"'-''-P,+»,]+SAB*SBA[pp*-pp*.pp*+pp*j 

(1+pp*)*<  1  Su'  *+l SBB 1  ,+l SAB 1  SBil  *> 


(1+pp  )span{[S(A,B)]} 


Th.«r»£or« 


Sp*a{[S'(A'.B')J>  .  ,pn{[S(i,B)])  .  imuiat 
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Appendix  E 

Derivation  of  tha  Mueller  Matrix f Ml  from  the  Scattering  Matrix  fSI 


Given 

h 


CE-1) 


and  so  that  : 

»0  *  1  *  IV*  +  ,hBl* 

*1  *  <3  "  I**1*  *  * 

tj  -  0  -  2B«<hih2*} 

tj  -  V  -  -Umlh^} 

satisfying 

*0*  "  *l1+s2a+g3*  m  z*  m  Q*-HJ*+Va 

tha  scattered  and  tha  incidant  Stokes  vactors  ga  and  g"  ara 
related  to  each  other  by  tha  Mueller  Matrix  [M]  : 

ga(A,B)  «  (M]g^(A,B) 

We  have  also  : 

ha  -  [S]  h1  ,  or 


V 

’SAA 

SAB 

1 - 

L 

m 

!«•  « 

SBA 

SBB 

m 

- 1 

_ 1 

Therefore  : 

'  ^A1  +  SAbV 

V  *  SW  *  W 


CE-2) 


CE-3) 


CE-4) 


CE-5) 


CE-6) 
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1.  Calculation  of  g-1,  g®  : 


iV  *  W* 


(SMll*i+SAflhBi)  '‘uWO 


S4ASAA*lhii|1+SASSAs*lhBt|‘+SUSAB\Si'’+SA3SAAV\l 

I ! SAA 1 1 '  I V  l: ‘+ 1  SAS  I  *  I », 1 1  *+SAASA*\  V‘ 

*  L  i*  * 

+C®U8AB  hAV  5 


Tharafora, 


lSM|,|h4l|‘*|SA»l‘|hB1|,*JR*{SMSAB\Sl*>  tW> 


But  wa  hava  : 


SAASAb\V*  ■  [*«<S4ASAB*,*JI“<S4ASAa*)lla*<1,AillBi*,+jI“{hA\i*)1 

*  H*<S4ASAB*)R.(ll41hJ1*)-X«{SAAS4j*}I.{hA\i*) 

+HR*{sMs4a*)ia{hA1iiB  *W.{hA1iji*}i»{s4AsAB*n 


Tharafora, 


®*<SAASAa\V*}  *  2S*<SAASAB*,R,<hAihB1*)'2Im<SAASAB*)I“{iA\i*t 


CE-8) 


Substituta  from  (E-8)  into  (E-7),  tharafora  : 


iV1*  "  !SAAl'|11A1|,+|SAS1,|hB1|,+2E,{SAASAB*,R*<hA\ii) 
-2Xo.{SaASaj*>Ib(1ia\1*) 


iVi*  ■  w* 


‘VaSbO^baVXbV^ 


(E-9) 
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'  SBa'  ’  I O  ,+l  SBB I  *  I  hBi>  ,+SBASBB  W^^BaSbXS1*!* 
I V  *  I  h4tl  1+l  SBB I  * 1 V I  ,+28*<SBASBB  W*1  (E'l0) 


Similiar  to  Equation  (E-8),  w*  can  writ*  : 


-  22«<SBASBB‘}R.{l.1\1*}-2Iin{SBiSBB*H.{hi\i*) 


(E- 11) 


Substituta  from  (E-ll)  into  (E-10),  than  : 

“  |SSAl»|hAi]^|SBBl^|hB1i'+;>K.{SBASBB*}R.{hA\i*} 

-21mi  SbaEBB*^ 

From  Eqs.(E-2)3  (E-9)  and  (E-12),  w«  can  writ*  g*  and 


(E— 12) 


g,  as  follows  : 

*g-  |h/|«*  lh,V 

-(|3AAl‘+lS3Al*)lhAl|l-K|SABIHSBBl‘)II>,1lI 

+2S*«SAASAB*+SBASBB*,R*{llAS1*>-2I"{SAASAB*+SBASBB*,I,"{hAi,1B' 


us*  th*  identity  : 

ax+by  *  i(a+b) (x+y)  +  i(a-b)(x-y) 


CE-13) 


CE-14) 


Tharafor*. 


*0S  ■  «lSAA|i+lSBA1‘+|SAB1,+|SBBi‘>ti,lAi|I*|llB1|’} 


*i{|SAAl,+ISBAl'-|SABl»-|SBBl*K|hA 

+2R.{S4ASAB*+SBASBB*)R.{bA\i*) 

-2I“<SAASAB%SBASBB*n“{hAV*> 


CE-1S) 
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Us  in#  thm  definition  (E-2)  into  (E-IS)  :  Th«ir®fors, 

*0S  *  “u<01+,,12*ll+,,13*21+!,14*3i 
wh*t'»  : 

-u-«|suI»+|sbaI1+|sm|*+|smI*) 

"12  ■«lSJ*+!SU|1-|S4sl''lS1)8l*} 

“l3  ■  R*<5MSAa*+SBASBS*) 


“24  *  WSUSAa%8BASBB*> 


CE-16) 


(1-17) 


4140  from  Eqs. (E-2, 9, 12)  ; 
gls  -  |liAa|a  -  \K*\* 

-  (lful*-|sBAl*),h4i|>+(|sia|*-|sBsl«)|hBi|* 

+iR“(SA4SAii*-SBASBB'')!(,<h4\i*)-H“{3AAS4B*-aB43BB*)I“{hA\1*J 


(E-liJ) 


0*ia#  th«  idmatity  (E-14)  with  (E-18),  ch«a  : 

*1*  ■  «lS4A!,-'SB4li+|S4B|,-|3BBl,>ClhAi|,+lhBi|,) 

*i(|SMl>-|SBAl*-|S4Sl**|SgBl*K|hA1|«-|ll31|*) 

'l'2S*^SAASAa  "SBASBB  * 

-2I“<SAASAB**aBASBB*>I”{hASi’’> 

*  “2At+“22*li+“23*2±+“24*3l 


where 


“21  -«ISAaI,-ISBA|,+,SABI*-|SBb',> 
a22  *  ^  S4iJ  SBA^  SAB^  1+ISB8^ 


(E-19) 


-  • .. M.'  I il':.1-;,. , I i . LV. Blj. All : [ .,1, 1 *1. J,' , i. li , .  JL.'l 


■M  *wr*i  m. 
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“23  *  to{SAASAB**SUSn*> 


“24  “  I“<SUSAa**SBASB8*) 


2.  Calculation  of  g.®  and  g3fl  : 


>aV*  “  <8uV+8akV)<sbaVHbV,') 

*»  ^  I  L.  ill.l.rt  lit  *1  L  n  ^fc.  it. 


SMSM  V^aAb  ^  « ^SAASBB 


hA\1  +SABi,BA  hA 


K  IK  i 


Th«rn£ct*o, 

»AV*  “  IR«tSMSBA’'1+'lI“i:S4ASBA*,l|hAil” 


4(B«<S43SM"H-JIn(SMaM*)]!UBi|‘ 

+(R«{Su3BB*)4jm{3MSBB*> )  [Sa{hA1hTt*)+i  Ia{h41hB1*}  1 
+(Ha{SAB3BA*)+jI«{3ABSBA*)HRa(h41hB1*)-jIo{liAiliBi*)]  (E-20 


From  (E-2)  and  (E-20)  : 
gj*  “  ESaUv^Uj**) 


"lA^BA  I*+Rb^SA3SBB  ^  * 


4Ra{hA^Ug  )S»{SAASBB  +3abSba  'SABSBA  ^ 

-  R«!SUSBA*+SAB31B*H  Ih^l  *+|hB1|  *> 


+R,(SUSBA"-SABSBB*KlhAil'-illBl|1> 
42Ra{SAASBB*+SABSBA*)Ro(hA^lB^*)-2Io(SAASBB*-SABSBA*}Iffl{liA^1B^*) 

'  “31*01+”32*lt-h"33*2l+”34*3i 


“31  '  a,{SAAS8A*+SABSBB*) 


wh.«r« 
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^  pi  ■mttrjv V  y-r  ,’fiT.t ' r«-r5i :*T.-f ^ 1  PTT  ,’,rJTwV 

wVi- 


*32  -  «.aAi3BA‘-3AJ3B!*> 


*33  “  S*£SM398*+SA33BA*) 


M-i  3,1  Lj{3 . . S--*"S . „S_ *} 
34  AA  BB  AB  BA  1 


Also  froa  (£>-2)  and  (£-20)  : 

a* 


(E-21) 


*3#  -  1 

-  -3(I«{3uSBA*Hhi1|*+I«{3AB3BB*)|hBt|‘ 

«.{h4\i‘'}I«(Su3BB*+SABSBA*)+Ia{hA\i*)R.{Su3BB*-3ABSBA*>] 

*  -I»t3AA3BA*+SAaSBa*K|llA1|>+|hBt|M 
-.uaC3AA3aA  l‘-|9j  I2) 

-2W31ASM*+8AB3BA*)R.{hAthB1*)-ai.(3AA3BB*-34B3JA*)Ia<llA\i*) 


whara  ; 


“41 


“42 


“43 


“44 


»  *  *4*1  „  „  i 

^41*0  +a42*l  *®43*2  *®44*3 


•I“<34i3BA*+84B3BB*1 


•l"<3aa3ba*-3absbb‘> 


-Ib(3u3bb*+8ab3bi‘1 

»»<3U3BB*-34B38A*> 


CE-22) 


la  summary,  than  tha  altmants  of  tha  Muallar  matrix  can 

ba  writ tan  in  tarma  of  tba  • lam ant 3  of  tha  acattaring  matrix  u 
follows  : 

(a)  Blatatic  caaa  : 

■u-i{|SAAl2+|SBAl2+lsABl2+|5BBl2) 

.  .  .  4,-  _  •  ■  "•'  11  Lb  '■  — ■ 
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Bi2-K!saaI1+ISbA,2-|sAbIj-,sbbI1> 


b13  *  Ra{SAASAB  +SBASBB  } 


b14  *  Ia<SAASAB*+SBASBB*} 

“>2i  *  ^  I SAA  ^  ■  SBA  ^  ^  SAB  ^  ^  SBB  ^  ^ 
«22-Ki^AAla*lSBA|,'1SAB,,+1SBBl#> 


n23  *  R**SAASAB  *SBASBB  * 
°24  *  Im*SAASAB  *SBASBB  * 
®31  *  R*^SAASBA  +SABSBB  } 
“32  “  R*^3AASBA  'SABSBB  } 
m33  *  R**SAASBB  +SABSBA  } 
®34  “  Ia*SAASBB  "SABSBA  * 
m41  “  ‘Im*SAASBA  +SABSBB  * 


n42  *  ’Ia<SAASBA  'SABSBB  } 


“43  *  -Io{SAASBB*+SABSBA*} 


°44  *  R*^SAASBB  ’SABSBA  * 


(E-23) 


b)  Monoatatic  cas*  : 

°U*H|SAA|i+2|SABll  +  |SBB|2> 
®12  *  ^^SAA^  *'*'SBB^^ 

®13  “  Rt^SAASA3  +SABSBB  * 


m14  *  Ia*SAASAB  +SABSBB  * 
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Si  *  -  a, 


21  av'“AA‘  '  BB1  4  “12 

n22  *  ^^SAA^l+'SBB^l"2^SAB^  ** 


®23  *  R*^SAASAB  ’SABSBB  * 


®24  “  ^SAASAB  ”SAESBB  * 


®31  *  R*{SAASAB*+SABSBB#>  “  “13 


“32  “  R**SAASAB  'SABSBB  *  *  ®23 


“33  “  Ra{SAASBB  >+lS4JBfl 


a34  »  Ia{S.  } 


BB 


“41 


,La*SAASAB  +SABSBB  *  *  "®14 


a 


42 


'WSUSil  "SASSBB  1  "  “°24 


"43  “  CSAASBB  1  '  '“34 


"44  *  R,£SUSBB  >-IS4bI‘ 
v«  also  do  notica  that 


“44  *  “33  +  °2?  '  “11 


G 


'■ . 
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Appendix  F 

Derivation  of  the  Modified  Mueller  Matrix  £Mm]  from  the 
Scattering  MatrixfSI  ; 

Dafina  the  modified  Stokes  parameters  for  the  polarization  vector 


A  A 

S  ■  hAhA  +  Vs 

cr-i) 

«*3  -  i*a'* 

-a 

H 

N 

V-'V 

Ami  *  2S*{h^hj  > 

■  u 

CF-2) 

®b3  ’ 

-  V 

and 


^(A.B)  -  [MmJs^CA.B) 
.We  liave  : 


V 

m 

•m  " 

SAA  SAB 

J 

a 

m 

SBA  SBB 

h1 

L  j 

m  « 

wm  — 1 

Therefore , 

"A*  *  +  *Alh»t 

V  ■  SBAhAi  +  33b41 

Then 


CP-3) 


CF-4) 


CP-5) 


'a5  “  l«*A*l*  -  »aV* 

*  fSAAiAi+SAah3t)CS4A\i*+SAs\i’) 

‘  'saaI  Vi^abI'I Vi'AAaVVAaSajYV* 


Therefore, 


V -iV* 


I SMI  *  |h/ 1  H  h,1!  ■+(sMsiS\\i*)+(sAAsAB*iiA\1*)' 

I  sul  *  I  h/ 1  ‘>|  su  |  *  |  h,1 1  *«s.<8usM \V*>  CT-4) 


But  we  have  : 


SAASAbW*  *  ‘R*{S4AS4a*HJIll{SAASAB*lltR*{hASi*)+JIo{h4\i*)1 

*  [R*{SAASAS*}R,{llAih3i*)-I“{3AASAB*)I"'(UASi*>! 

+j  [  WSMSAa*)8«(hAillBi*}+R.{SAAS4s*}Im{liA\i*)  1 


Thar afore, 


WyAS1*1  *  2S*{3AASAa*)R“{hA\1*)‘2I“{3AASAB*)lB{hA\1*) 


or-?) 


^VAS1*1  *  2I“<SAASAS*JRe(hA\i*)+2S*{SAASA3*)In{hA\i*1 


CF-8) 


Subatituta  from  (F-7)  into  (F-6),  therefore  : 

V-  'Vi* 


■  I  SiAl 1 1  hAi|  »+|  Sa3  H  hj,1 1  >+2S.{SAASA3*)R.{hA\i*} 
-2Ln{SAASAJ*}Ia{hA\i*) 


Therefore, 


:AS  *  lSul,IAi+|SA3|1I8i+R*{SAASAB'*,Di+I"<SAASAB*,vi 


also  wa  have  : 


V“  'V*  *  w* 


<SBAhA  +SBBh3i)CSBA  hA  +SBB  hBi  5 
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■  I V  '  I  “a1!  ,+l  SBB I 

■  I SBAI 1 1  u/l  *+|  SBB|  »  I  hj1!  *+2Sa{SBASBB*)R.{hA1hBi*> 
-2Ia{SBASaB*)Ia{hA^iBi*} 


Than 

V  ’  lSBAl,IAi+lSBBl,IB1+R,tSBASBB*>0i+I"<SBASBB*>vi 

Also  wa  hava  : 

\  L\z  \  **> 


(II) 


“aV*  *  <SAARaHbV><SBaV"+SBbV~> 


*  'uOO  ^3AaSBa*lbBi|  ,+SAASBb\  V*+3ABSBA*<hAY*3 

*  R*<  3AASBA*>  I  ''A1 1  J+R*<  SABSB8*>  I ' O  * 

«*JI«(SabS,b*>  Ih,1!  * 

+(S»(SAASBB*>+JI“{SAASBB*>liR*{hA\1*}+JI“<hA\i*» 

*IE*<SABS3A*,+JIm<SABSBA*lltR*<hA\1*>'-lIm<hAlh31*)I 


Than 


a,.  3*  a 


2S.{1>a  )  -  2S.{SaaSba  )lhAl*+2R.{SABSBB  >11131* 

«Re{hAy*>R.<S4ASBB\SiBSBA*> 

-2Ia(hA\i*>In,(5AASB8*.SAB3BA*> 

Tharafora, 

O3  -  2B«(SAASBA*>IA1+2R.{S1DS„*>I„i 


AB  BB  '  B 


*R.{SAASBB*+SABSBA*HJi+Ia(S4ASBB*-SABSBA*)Vi 


(III) 


■  I 


HttdiaiiMitiaiitaiiMtUaMiitfiiiMaaaMiMi 


-  -2I»{S4iSM*)|hii|‘-2I«{SASSBS*)|hBl|* 


t 

i 


f 


-2S.{SA4SM*-SiB3BA*)L.<hA1ll3i*' 

Tharafora, 

t 

'WSAA3BB*+Sia3Bi*lai+R,<34ASBa‘-34aSBi*,V'1 


la  surasary,  from  (I)  to  (IV),  wa  can  write  [Mm]  an  : 
r 


l<W* 

r*^sAASAB  j 

W*aAj*> 

|sbbI* 

S*{SBASBB*> 

* 

»*«u*u*> 

*®*^sabsbb  * 

r*{saasbb  kSabsba  } 

XaiS^Sjg  'S^3S3A  > 

-2«S4BSBB*} 

~-1®' °MSBB  +SABSBA  } 

r*^saasbb  "sabsba  * 

(V) 


1 

) 
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Appudix  G 

Derivation  of  the  Scattering  Matrix  f  3 1  from  the  MusHer 
Matrix f Ml  : 

It  is  required  to  calculata  the  scattering  matrix  [S]  from  the 
Mueller  matrix[M]  in  the  general  bistatic  and  sonostatic  cases. 

The  scattering  matrix  [S]  relates  the  scattering  and  the 
incident  polarization  vectors  by  the  equation  : 

hS(A,B)  -  [S]  hi(A,B) 
where  : 


m 

sab 

m 

"'3aa' 

ej*AA 

1 - 

3 

♦ 

•n 

e 

to 

2 

w 

_ i 

SBB 

'SBA 

e^BA 

sBB!aJ*B8_ 

Also  the  Mueller  matrix  relates  the  scattered  and  the  incident 
.Stokes  vectors  by  : 

ga(A,B)  -  [Ml^CA.B) 
where  : 


*0 

&  * 

*1 

* 

*2 

OaO^h,*} 

*3 

•  m 

-21^(^112*) 

1.  Calculation  of  the  magnitudes  I  S^| ,  |  S^j ,  |  Sj^l  and  |SBBI  : 


From  Appendix  (E) ,  we  have  : 

«U-HlSAAiHSBAl^sABIHSBBla> 


(8-1) 

(8-2) 

(8-3) 

(8-4) 


(8-5) 


I 


“l2  "  ^  lSAA^  1"‘S4S^ 

Vi(|34Al*-|SBBl‘+|SABl,-lSBBl,) 

aj2-«|Sj»-|SB4l*-|SA1>l**|SMl») 

n  h*v«  : 


*11  +  *12  *  |SU1’  +  ^BA1* 


*21  +  *22  *  I®**1*  ’  ^BA1* 


than 


and 


*11  +  *12  +  *21  +  *22  "  2|V“ 


*11  +  *12  ‘  *21  '  *22  ’  J|SBAI*  ’  Cll*r,f0r* 
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CG-6) 

(G-7) 

(G-8) 

(G-9) 

(G-10) 


ISaaI  * 

iC“ll4*‘l24*214,"22) 

(G-Xl) 

|sba'  “ 

i(*u4»i2-«2i-«22) 

CG-12) 

Siailiary  : 

BU  *  “12 

-  I®a»I*  +  |3bbI* 

“21  *  “22 

■  ISab1’  •  l®nl* 

Thor • for a. 

!sab!  ’ 

i(“ir*12+*2l‘"22) 

(G-13) 

Isbb' 

i(*U-*12'*21+*22) 

(G-14) 

2.  Calculation*  of  tha  phasos  $BB  : 

H 


Also,  wo  havo  from  Appandix  (E)  : 


*13  ”  S“(SAASAB  *  SBASCS  J 


*14  "  I“<3AA3AB*  *  SBASBB*) 


i- '  J  i: 


t  >„  --  •BJT-.r_wnfii-i  ‘  tl1  1  '  r  - - rasa*?  -  TJL..E.  li-.'i-VJT 
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13m  : 

■a  *  *“».  “  yu'  *  *u*M* 

ilM, 

■a  ■  wsu*ii  - 

•»  *  «Vu*  *  W> 

•a  *■  i"»  *  Vu'  *  SjaS** 

fra  CS-U)  ad  CS-14)  : 

<*u*i»x4>  *  C«23+1^A>  “  ausa* 


<*U**B3  *  K*w*«ifc>  *  2-1 3^*1 
Stmfaw, 

■  "14*^4. 

•u  *  «u+«aaBC— — ) 

*u**a 

Um  tna  (S-m  ad  (s-u)  . 

‘"u*^  *  <«a*J-z*.>  •  W 
c«u-«a>  *JCw)'»2l»SAll^,Ult*a4'*u0 

*14-*2V 

*u"*u 

Tram  ifpmdix  £,  m  41m  hm  : 

•n  *  ^suSA#**iaV> 


(5-13) 


(5-16) 


(5-17) 


CG-li) 
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m41  =  ■Im{SAASBA  +SABSBB  } 


°32  =  Re{SAASBA  'SABSBB  } 


m42  *  "Im{SAASBA  'SABSBB  } 


Therefore, 


*  * 

m3l"jm41  =  SAASBA  +SABSBB 


and 


=  8aa8ha*-sabsbb*  '  then 


j2  J  42  AA  BA  AB  BB 
(ra3rJm41)'‘<-m32'jra42)  =  2SABSBB 


* 


(ra3rm32)+j(m42'”4i:)  =  2|SAB1|SBB|eK0AB  ^ 


CG-19) 


(G-20) 


therefore 


“42'”41 

*ab"*bb  ■  arotan  ( - : 

m3l"m32 

m4l"°42 

*BB  =  *AB+arctan(: - } 

m3l‘ra32 

From  (G-18),  (G-21) ,  therefore  : 

m14‘m24 

*BA  =  ^B+arCtan( - 5 


ml3"ai23 
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ffi14'm24 


m4rm42 


*BA  =  ^AB+arCtan(' 


-)+arctan( - ) 


m13_m23 


m3l‘m32 


(G-22) 


In  the  monostatic  case  : 
From  Appendix  (E),  we  have  : 

m21  =  ”12 

m31  =  “l3 

ra32  =  m23 

m41  =  "ml4 

m42  =  “m24 
m43  =  -m34 


m44  =  m33+ra22_m^2. 

Therefore, 

^AA1  =  i(mll+2m12+ta22) 

^SAB^  =  ^(mll'm22)  =  *SBA^ 

I^B1  =  *(Bir2B12+B22) 

ml4+°24 


*AA  =  *AB+arctan(- 


0ab  is  arb. 


m14'm24 


-m14+°24 


*SA  =  ^AB+arctan(' 


-)+arctan(- 


m13’m23 


m13*m23 


AB 


-7F  ."W  -i  A  mWJUHJHMJI  IM^L 
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m14‘°24 

*BB  =  ^AB'arCtan( - } 

m13‘m23 

In  suromary  : 

1.  Bistatic  case  •, 

lSAA*  iC“,ll+’,12+“21+B22) 

ISAB1  * 

lSBA^  =  \J  I  (mll+m12"“2l'“22'1 

ISflB1  *\/40"  n'*X2'l‘zi*m22> 

”l4™24 

♦M  =  *AB+atCtan( - J 

$  is  arbitrary  (may  be  =  0) 

AB 

ml4'm24  m4r°42 

0  =  0  +arctan( - )-farctan(  *  ) 

»13-"23  -«-« 


m4l'm42 

*BB  =  *AB+arCtan( 

m3l"*32 


2.  Monostatic  case  : 

^AA1  =,\/r*(“‘ll+2B124l"22) 

^AB1  ‘V  iC“u'“22)  =  1 SBA^ 
ISbb1  -J  i(oll‘2B12'h,22) 
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m14+ffi24 

*AA  =  *AB+arCt£m( - } 

ro13+m23 

*ab  is  arb* 

0BA  =  *AB 

®14““24 

*BB  =  *AB'arCtan( - > 

m13‘m23 


Appendix  H 

Derivation  of  the  Scattering  MatrixfS]  from  the  Modified  Mueller 
MatrixfMml  : 


From  Appendix  F,  Eq.  (V),  we  have 

Mu  •  |smI‘ 

M12  -  'SAbI‘ 


H13  *  Ra(SAASAB  > 


M14  “  ImtSAASAB  1 


M21  =  SBA^ 


M22  =  ISBb'* 


«23  “  R«(SBASBB  > 


M24  ■  «SBASBS  > 


M3l  *  2E.(SmSba  } 


M32  -  2Re(SABSBB  > 


M33  Ra(SAAS8B  +SABSBA  ! 


”34  Ia{SAASBB  'SABSBA  1 


“42  -  -2«SAASBA  > 


"42  ■  -2I"<SABSBB  > 


M43  “  'IlfSAASBB  +SABSBA  1 


M44  =  RetSAASBB  "SABSBA  1 


Calculation  of  the  amplitltudea  : 
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ISAb‘  =\ph2 

>sba'  -^T 

^SB»'  =  \j  M22 

Calculation  of  the  phases  : 

Mi3+JMu  = 

”  SAASAB 

=  !SAAllSAB|eJ<*M'*AB) 

Therefore, 


M 


14 


*“WAA'W  = 


M 


13 


Then, 


M 


14 


*AA  "  ^AB+arctan(' 


M 


13 


also, 


M34+M43  “  Im{SAASBB  "SABSBA  J"Im(SAA3BB  +SA8SBA  } 

-  ImCS^Sgg  > - Im< SABSBA  }-ImiSAASBB  }'Im{3ABSBA  } 


=  -2IoKSabSba  > 


M33'M44  RetSAASBB  +SABSBA  *‘Re’SAASB3  "SABSBA  * 


*  2Re{SABSBA  > 


and 
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then 


(M33'M44)"j(M34+M43) 


2 [^e{SARSpA*}+jIm{SABSBA*} 


=  2sabsba 


AB  BA 
* 


-  2|SAJl|SBAleJ(*AB-*BA) 


M34+M43 


tan(^AB-^gA)  - 


Therefore, 


M33‘M44 


M34+M43 


*AB“*BA  =  -arctan(- 


M33'M44 


or 


M34+M43 


*BA  =  ^AB+arCtan  (' 


”) 


M33"M44 


also  we  have 
*32*^2 


M32‘-iM42  "  2tRe<SARST,/>+jIra{SABSBB*^ 


AB  BB 

=  2sabsbb 


=  2lS4BIISBB|eJUAB">SB) 


then 


M 


42 


^♦ab'W  =  • 


M 


32 


or 


M42 

^BB  =  *AB+arctan( - > 

M32 
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In  summary  :  for  the  bistatic  case  : 


and 


*14 

»AA  =  *AB+arctan( - 5 


M 


13 


<^2  is  arbitrary  (may  be=0) 


*34+*43 


>BA  =  ♦AB+«ctan(- 


*33~M44 


*BB  =  «AB+arCtan( - 5 

*32 

simplify  in  the  monostatic  case  so  that  IS^MS^I  and 
*AB^*BA’ 


APPENDIX  (I) 
COMPUTER  PROGRAMS 
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This  Appendix  includes  three  computer  programs .  The  first  one 
calculates  the  Mueller  matrices  [M]  and  [Mm]  from  the  Scattering 
matrix  [S] .  Also,  it  calculates  the  C0P0L  and  XPOL  nulls  and 
their  representation  on  the  Poincare  Sphere. 

The  second  program  reconstructs  [S]  with  relative  phase  from 
[M]  or  [Mm] . 

The  third  one  reconstructs  also  [S]  from  the  knowledge  of  the 
spherical  coordinates  on  the  Poincare  Sphere  of  either  two  C0P0L 
nulls  or  one  COPOL  and  one  XPOL  null  . 
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1-1  FIRST  COMPUTER  PROGRAM  : 


This  program  calculates  the  Mueller  matrices  [M]  and  [Mm]  f  ->ra  the  scattering 

matrix  [S]  in  the  bistatic  and  monostatic  cases  using  the  Equations  of 

Section  3.6.1.  Also,  it  calculates  the  COPOL  and  XPOL  nulls  and  their  representation 

on  the  Poincare  sphere  using  the  equations  of  Section  3.5. 

//BAKRY  JOB 
/*JOBPARM  R=348 
//  EXEC  WATFIV 
//SYSIN  DD  * 

$ JOB  NOEXT 
C 

C  THIS  PROGRAM  CALCULATES  THE  MUELLER  AND  THE  MODIFIED  MUELLER 
C  MATRICES  FROM  '.HE  SCATTERING  MATRIX  IN  BISTATIC  AND  MONOSTATIC 
C  CASES. 

C  ALSO,  IT  CALCULATES  THE  COPOL  AND  XPOL  NULLS  FOR  S  AND  REPRESENT 
C  THEM  ON  POINCARE'  SPHERE  . 

C  PUT  LANE=0  IF  SIGMA'S  AND  PHI'S  ARE  KNOWN 
C  PUT  LANE.NE.O  IF  S(I,J)  IS  KNOWN 

C 

COMPLEX  S(2,2),RH0(2),CI 

DIMENSION  AM(4, 4) ,AMM(4,4) , THETA (2),PHI(2) 

COMMON  Cl, PI, PI 1, PI 180 
CI-CO. 0,1.0) 

PI=4.0*ATAN(1.0) 

PI1=180.0/PI 

PI180=PI/180 

C 

C  DATA 

C 

NDATA=9 

LANE=0 

DO  999  11=1 ,NDATA 
IF(LANE.EQ.O)  GO  TO  160 
READ  55,((S(I,J) , J=1 ,2) ,1=1,2) 

55  F0RMAT(4E10.3) 

GO  TO  161 

160  CONTINUE 

C  SIGHH,  SIGHV,  SIGW  IN  DECIBELS 

C  PHIHH.PHIW  IN  DEGREES 

C 

READ  (5 , 155  )  ASPECT ,  SIGHH ,  SIGHV ,  SIGW ,  PHIHH ,  PHIW 

155  F0RMAT(6F10.2) 

WRITE(6, 156)  ASPECT,  SIGHH,  SIGHV, SIGW,  PHIHH,  PHIW 

156  FORMATC'l' ,' ASPECT  ANGLE  =  ' ,F1G . 2 , 2X, 'DEGREES '//' O' ,' SIGHH  =  ' 

1 ,F10. 2,2X, 'DECIBELS' / ' O' , ' SIGHV  =  ' ,F10 . 2 , 2X, ' DECIBELS '/' 0 ' , 

2 ' SIGW  =  ’ ,F10.2,2X, 'DECIBELS '/' O' , 'PHIHH  =  ' ,F10 . 2 , 2X, ' DEGREES ' 

3/  '  0 '  , '  PHIW  =  '  ,F10. 2, 2X, 'DEGREES'//) 

CALL  SCAT ( SIGHH ,  SIGHV ,  S IGW , PHIHH ,  PHI W ,  S ) 

161  CONTINUE 


J 


ooo  o  ooo  o  n  o  o  o 
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S  C  2 , 1)=S(1,2) 

PWR=NORMALI ZED  FACTOR 


PWR=CABS(S(1,1))**2-K:ABS(S(2,2))**2+CABS(S(1,2))**2+CABS(S(2,1))** 
12 

PVT2=SQRT(PWR) 

WRITE(6 ,50) 

50  FORMAT (' 1 THE  SCATTERING  MATRIX  IS  :'//) 

DO  1  1*1,2 

WRITE(6 ,60) (S(I,J) ,J=1,2) 

1  CONTINUE 

60  FORMATC'O' ,2X,2(2E10. 3,5XP 

WRITE  NORMALIZED  SCAT).  ,‘R  J  MATRIX 
WFITE(A ,51) 

51  FORMATC'O' , 'THE  NORMALIZED  SCATTER T.NG  MATRIX  IS  :'//) 

DO  11  1=1,2 

: 1  WRI'  ~ (6 , 60) (S(I , J)/PWR2 , J=1 ,2) 

CALL  STOM(S.AM) 

WRITE(6 ,70) 

*/0  FORMAT  ('•■', 'THE  MUELLER  MATRIX  IS  :’//) 

DO  2  1=1,4 

WRITE (6, 80) (AM (I , J) ,J=1,4) 

2  CONTINUE 

80  FORMATC'O' ,2X,4(E10 . 3 ,5X) ) 

XX=AM  (3,3)  -irAM  (2,2)-AM(i,l) 

WRITE ( 6 , * ) XX , PVR , PWR2 

WRITE  NORMALIZED  MUELLER  MATRIX 

WRITE(6 , 71) 

7 i.  FORMATC'O' , 'THE  NORMALIZED  MUELLER  MATRIX  IS  :'//) 

DO  12  1=1,4 

12  WRITE (6 ,80) (AM(J , J)/PWR, J=l,4) 

CALL  STOMM(SrAMM) 

'WRITE  (6, 90) 

90  FORMATC'l' ,'THE  MCIIFIED  MUELLER  MATRIX  I'S  :'//) 

DO  3  1=1,4 

WRITE (6 ,80) (AMM(T , T) , J=l,4) 

3  CONTINUE 
C 

C  WRITE  NORMALIZED  MODIFIED  MUELLER  MATRIX 
C 

WRITE (6 ,91) 

91  FORMATC'O' ,'THE  NORMALIZED  MODIFIED  MUELLER  MATRIX  IS  :'//) 

DO  13  1*1,4 

1 3  W R 1 IK ( 6 , 0  0 ) ( AMM ( I , J ) / PWR , J= 1 , 4 ) 

C 

I.SW*0 

CALI.  COPOLCS  ,RH0, THETA, PHI , ISW) 
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IF(ISW.EQ.l)  GO  TO  120 
WRITE (6, 100) 

100  FORMATC'l' , 'THE  COPOL  NULLS  ARE  :’//) 

DO  4  1=1,2 

WRITE (6, lQl)RHO(I) , THETA (I) ,PHI(I) 

4  CONTINUE 

101  FORM^C’O' ,  'RHO  =  '  ,2G14. 7// ' 0 '  ,  'THETA  =  '  ,F8. 3, 2X,  'DEGREES'// 
1 ' 0 ' , *  PHI  =  ' ,F8. 3, 2X, 'DEGREES'//) 

120  CONTINUE 

CALL  XP0L(S, RHO, THETA, PHI, ISW) 

IT'(ISW.EQ.2)  GO  TO  999 
WRITE(6, 110) 

110  FORilATC'O'  ,'THE  XPOL  NULLS  ARE  :'//) 

DO  5  1=1,2 

WRITE (6 , lOl)RHO(I) ,THETA(I) ,PHI(I) 

5  CONTINUE 
999  CONTINUE 

STOP 

END 

C 

SUBROUTINE  STOM(S,AM) 

C  THIS  ROUTINE  CALCULATES  THE  MUELLER  MATRIX  AM  FROM  SCATTERING 
C  MATRIX  S  IN  BISTATIC  AND  MONOSTATIC  CASES. 

C  G=(I,Q,U,V) 

C 

COMPLEX  S (2 , 2) .CONJG , SI , S2 
DIMENSION  AM(4,4) 

SAA=CABS(S(1,1))**2 
SAB=CABS(S(1,2))**2 
SBA=CABS  (.S  (2 , 1)  )**2  ■ 

SBB=CABS(S(2,2))**2 

AM ( 1 , 1 ) = ( S AA+S AB+S B A+SB B ) / 2 . 0 

AM(l,2)=(SAA-SAB+SBA-SBB)/2.0 

AM(2 , 1)=(SAA+SAB-SBA-SBB)  /?. .  0 

AM(2,2)=(SAA-SAB-SBA+SBB)/2.0 

S1=S(1,1)*C0NJG(S(1,2)) 

S2=S (2 , l)*CONJG(S (2 , 2) ) 

AM(1 ,3)=REAL(S1+S2) 

AM(1,4)=AIMAG(S1+S2) 

AM(2 ,3)=fuiAL(Sl-S2) 

AM(2,4)=AIMAG(S1-S2) 

Sl=S(l,l)*CONJG(S(2,l)) 

S2=S(l,2)*CONJG(S(2,2)) 

AM(3,1)=REAL(S1+S2) 

AM(3,2)=REAL(S1-S2) 

AM(4,1)=-AIMAG(S1+S2) 

AM(4,2)=-AIMAG(S1-S2) 

S1=S(1,1)*C0NJG(S(2,2)) 

S2=S(1,2)*C0NJG(S(2,1)) 

AM(3,3)=REAL(S1+S2) 

AM(3,4)=AIMAG(S1-S2) 

AM(4, 3)=*AIMAG(S1+SP) 

AM(4,4)=REAL(S1-S2) 

RETURN 


-  v,\  al'fl  u* 


r  i 


---  i- 
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END 

c: 

SUBROUTINE  STOMM(S ,AMM) 

C  THIS  ROUTINE  CALCULATES  THE  MODIFIED  MUELLER  MATRIX  AMM  FROM  THE 
C.  SCATTERING  MATRIX  S  IN  BISTATIC  AND  MONOSTATIC  CASES. 

C  GM=(IA,IB,U,V) 

C 

COMPLEX  S(2,2))CONJG,Sl,S2 
DIMENSION  AMM (4, 4) 

AMM(l,l)aCABS(S(l,l))**2 
AMMC 1 , 2)=CABS (S ( 1 , 2) )**2 
S1=S(1 , 1)*C0NJG(S(1 ,2) ) 

AMM(1,3)=REAL(S1) 

AMM(1,4)=AIMAG(S1) 

AMM (2, 1)=CABS (S(2 , 1) )**2 
AMM  (2 , 2)=CABS  (S  (2 , 2) )**2 
S1=S(2, 1)*C0NJG(S(2 ,2) ) 

AMM(2,3)=REAL(S1) 

AMM(2,4)=AIMAG(S1) 

S ] =S ( 1 , 1 ) *C0N JG ( S (2 , 1 ) ) 

AMM(3,1)»2.0*REAL(S1) 

AMM(4,1)=-2.0*AIMAG(S1) 

S1*S(1,2)*C0NJG(S(2,2)) 

AMM (3 , 2)=2 . Q*REAL(S1) 

AMM(4,2)=-2.0*AIMAG(S1) 

S1=S(1, 1)*C0NJG(S(2,2) ) 

S2=S(1,2)*C0NJG(S(2,1)) 

AMM(3,3)=REAL(S1+S2) 

AMM(4,3)“-AIMAG(S1+S2) 

AMM(3 ,4)=AIMAG(S1-S2) 

AMM(4,4)=REAL(S1-S2) 

RETURN 

END 

C. 

SUBROUTINE  COPOL(S ,RHO , THETA, PHI , ISW) 

C  THIS  ROUTINE  CALCULATES  THE  COPOL  NULLS  AND  REPRESENT  THEM  ON  THE 
C  POINCARE'  SPHERE  . 

C 

COMPLEX  S (2 , 2) ,RHO (2) ,A,B ,C,D,CSQRT 
DIMENSION  THETA(2) ,PHI (2) 

A=S(2,2) 

B=S(1,2)+S(2,1) 

C-S (1,1) 

D=CSQRT(B*B-4.0*A*C) 

ERROR* l.OE- 10 

WRITE(6,*)A,B,C,D 

ERRB«CABS(B) 

ERR=CABS(A) 

IF (ERR. LE. ERROR. AND. ERBB.LE. ERROR)  GO  TO  30 
IF (ERR. LE. ERROR)  GO  TO  10 
RH0(1)*(-B+D)/(2.0*A) 

RH0(2)=(-B-D)/(2.0*A) 

GO  TO  20 

10  WRITE (6, 50) ERR 
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50  FORMATC'O’ , ’WE  HAVE  ONE  ROOT  BECAUSE  A  =  ',G14.7/) 
RHO(l)=-C/B 
RH0(2)=RH0(1) 

20  CONTINUE 

CALL  POINCRCRHO, THETA, PHI) 

GO  TO  1 
30  CONTINUE 
ISW=1 

WRITE(6 ,40) 

40  FORMAT ('O’ ,'A  =  B  =  0.0  ,  NO  ROOTS  . '/) 

1  CONTINUE 
PEYURN 
END 
C 

SUBROUTINE  XPOL ( S , RHO , THETA , PHI , ISW) 

C  THIS  ROUTINE  CALCULATES  XPOL  NULLS  AND  REPRESENT  THEM  ON  THE 
C  POINCARE’  SPHERE  : 

C 

COMPLEX  S ( 2 , 2 ) , RHO ( 2 ) ,A,B,C ,D ,CONJG,CSQRT 
DIMENSION  THETAC2) ,PHI(2) 

A=S(2 , 2)*CONJG(S (2, 1) )+S (2 , 1)*CQNJG(S (1,1)) 
B=S(2,2)*C0NJG(S(2,2))-S(1,1)*C0NJG(S(1,1)) 

B=B-S(1 ,2)*C0NJG(S(2 , 1) )+S (2 , 1)*C0NJG(S( 1 ,2)) 

C=S (1 , l)*CONJG(S (1 , 2) )+S (1 , 2)*C0NJG(S(2, 2) ) 

C=-C 

D=CSQRT(B*B-4.0*A*C) 

ERROR-1. 0E- 10 
ERRB=CABS(B) 

ERR=CABS(A) 

WRITE(6,*)A,B,C,D 

IF ( ERR. LE. ERROR. AND. ERRB.LE. ERROR)  GO  TO  30 
IF (ERR. LE. ERROR)  GO  TO  10 
RHO(1)-(B+D)/(2.0*A) 

RH0(2)=(B-D)/(2.0*A) 

GO  TO  20 

10  WRITE(6,50)ERR 

50  FORMATC'O' ,'WE  HAVE  ONE  ROOT  BECAUSE  A  =  ' ,G14 .7/) 

RH0(1)=C/B 

RH0(2)=RH0(1) 

20  CONTINUE 

CALL  POINCRCRHO, THETA, PHI) 

GO  TO  1 
30  ISW=2 

WRITE (6 ,40) 

40  FORMATC’O' ,'A  =  B  3  0.0  ,  NO  ROOTS  . '/) 

1  CONTINUE 
RETURN 
END 
C 

SUBROUTINE  POINCRCRHO, THETA, PHI) 

COMPLEX  RH0C2) ,CI ,Q,CONJG,YY 
DIMENSION  THETA(2) , PHI (2) 

COMMON  Cl, PI, FI 1, PI 180 
AMIN=1.E-10 
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DO  1  1=1,2 
YY=1.0+CI*RH0(I) 

YY1=CABS(YY) 

IF(YYl.LE.AHIN)  GO  TO  10 
Q=(1.0-CI*RH0(I))/YY 
Q2=CABS(Q)**2 
Q3=-(1.0-Q2)/ (1.0+Q2) 

W1=AIMAG(Q) 

W2=REAL(Q) 

IF(ABS(W1) .LE. AMIN. AND. ABS(W2) .LE.AMIN)  GO  TO  40 
PHI (I)=-ATAN2(W1 ,W2) 

PHI(I)=PHI(I)*PI1 
GO  TO  20 
40  WRITE (6 ,50) 

50  FORMATC'O' , ’PHI  IS  ARBITRARY  ,  REAL  AND  IMAG  Q  ARE  ZERO*/) 
GO  TO  20 
10  CONTINUE 
Q3-1.0 
WRITE(6 ,30) 

30  FORMATC'O’ , 'Q  IS  INFINITY  ,  PHI  IS  ARBITRARY  '/) 

20  TKETA(I)=ARC0S(Q3) 

THETA ( I ) =THETA ( I ) *P I 1 
1  CONTINUE 
RETURN 
END 

SUBROUTINE  SCAT(SIGHH,8IGHV,SIGW,PHIHH,PHIW,S) 

COMPLEX  Cl ,  S (2 ,2) 

COMMON  Cl ,PI ,PI1 ,PI180 
PRINT, 'Cl  ,  PI  =  ' ,CI,PI 
PHIHH=PHIHH*PI180 
PHI  W=PHIVV*PI  180 
SIGHH=SIGHH/20 . 0 
S1-10.0**SIGHH 

S ( 1 , 1)=S1*(C0S(PHIHH)+CI*SIN(PH1HH) ) 

SIGHV=SIGHV/20.0 
S2=*10 . 0**SIGHV 
S(1,2)=S2 
S(2,1)=S(1,2) 

SIGVV=SIGW/20 . 0 
S3=10.0**SIGW 

S(2,2)-S3*(C0S(PHIW)+CI*SIN(PHIW)) 


PRINT, 1 

RETURN 

END 

'  MAG (S ( 1 , 1)) 

,MAG(S 

(1,2)), MAG ( S ( 2 

,2))  * 

’  ,S1  ,S2,! 

$F.NTRY 

0.0 

3.1 

-24.4 

4.6 

203.5 

193.5 

9.50 

0.5 

-14.9 

2.1 

170.2 

-156.2 

26.0 

3.6 

-16.1 

2.8 

128.0 

121.0 

42.75 

1.2 

-29.5 

1.2 

65.1 

75.6 

77.71 

21.5 

-18.0 

21.5 

-8.3 

-10.8 

96.0 

-0.5 

-22.1 

-0.8 

-36.0 

-45.5 

108.0 

-23.2 

-24.8 

-5.1 

108.3 

-162.2 

125  50 

-1.4 

-28.5 

0.6 

-22.9 

-60.4 

172.0 

17.3 

-15.2 

18.4 

133.4 

-131.4 
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1-2  SECOND  COMPUTER  PROGRAM  : 

This  program  reconstrucs  the  scattering  matrix  [S]  with  relative  phase  from  [M] 

or  [Mm]  using  the  equations  of  Section  3.6.2. 

//BAKRY  JOB 
/*JOBPARM  R=348 
//  EXEC  WATFIV 
//SYSIN  DD  * 

$J0B  NOEXT 
C 

C  THIS  PROGRAM  CALCULATES  THE  RELATIVE  PHASE  SCATTERING  MATRIX  (S) 

C  FROM  MUELLER  (M)  AND/OR  MODIFIED  MUELLER  (MM)  MATRICES  IN  BISTATIC 
C  AND  MONOSTATIC  CASES. 

C 

C  DATA  :  MDATA  =  NUMBER  OF  CM) *  S  MATRICES  TO  BE  ENTERED  AS  DATA  . 

C  MMDATA  =  ,,  , , (MM) 'S  ,,  ,,  ,,  ,  ,  ,,  ,,  . 

C  M0N=0  FOR  MONOSTATIC  ,  MON  <>  0  FOR  B I STATIC 
C  NROW=  NUMBER  OF  ROWS 

C 

COMPLEX  S(2,2),CI 

DIMENSION  AM(4f4),AS(2,2) , PHI (2,2) 

COMMON  PI, PI 1, Cl, MON 
I?  I  *4 . 0*ATAN  (1.0) 

PI1-180 . 0/PI 
CI=(0. 0,1.0) 

C  WRITE(6,*)PI,CI 

M0N=0 
NR0W=4 

IF(MON.EQ.O)  NR0W=2 
MDATA  *  3 
MMDATA  =  3 

NDATA  a  MDATA  +  MMDATA 

DO  999  I 1=1, NDATA 

READ  55  ,  ( (AM(I ,  J) ,  J=1 ,4) 1=1  ,NR0W) 

55  FORMAT (4E 10. 3) 

IF (I I. GT. MDATA)  GO  TO  10 
WRITE(6,50) 

50  FORMAT ('1',' THE  INPUT  DATA  OF  THE  MUELLER  MATRIX  IS  : '//) 

DO  1  1*1 ,NRQW 

1  WRITE (6, 60) (AM(I , J) , J*1 ,4) 

60  FORMAT ('O' ,2X,4(E10. 3,5X)) 

CALL  AMTOS(AM,AS ,PHI) 

GO  TO  20 
10  WRITE(6 , 70) 

70  FORMAT (' l’ , 'THE  INPUT  DATA  OF  THE  MODIFIED  MUELLER  MATRIX  IS  :'//) 

DO  2  1=1 ,NROW 

2  WRITE (6, 60) (AM (I , J) , J=1 ,4) 

CALL  AMMTOS(AM,AS,PHI) 

20  CONTINUE 
DO  30  J»l,2 
DO  30  1=1,2 

SCI ,  J)=AS(I ,  J)*(COS(PHI  (I ,  J)  )+CIlVSIN(PHI(I ,  J) ) ) 
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PHI (I , J)=PHI (I , J)*PI1 
30  CONTINUE 

WRITE (6,80) PHI (1,2) 

80  FORMAT(' , 'THE  SCATTERING  MATRIX  WITH  ARBITRARY  PHASE  PHIAB  =  * , 
1F8.3.2X, 'DEGREES'//) 

DO  3  1=1,2 

3  WRITE (6, 60) (S(I,J),J=1,2) 

WRITE (6, 90) 

90  FORMATC'-' , 'THE  SCATTERING  MATRIX  IN  POLAR  FORM  IS  : '//) 

DO  4  1=1,2 

4  WRITE (6, 60) (AS (I , J) ,PHI (I , J) ,J=1,2) 

999  CONTINUE 

STOP 

END 

C 

SUBROUTINE  AMTOS(AM,AS,PHI) 

C 

C  THIS  ROUTINE  CALCUUTES  THE  RELATIVE  PHASE  SCATTERING  MATRIX  FROM 
C  MUELLER  MATRIX  IN  BISTATIC  AND  MONOSTATIC  CASES  . 

C  PHI (1,2)  IS  ARBITRARY  E.G.  (  =  0.0  ) 

C 

COMPLEX  Cl 

DIMENSION  AM(4,4),AS(2,2) , PHI (2, 2) 

COMMON  PI, PI 1, Cl, MON 
C‘  WRITE(6  ,*)PI  ,CI 

C. 

C  CALCULATE  THE  MAGNITUDES  OF  (S) 

C 

IF(MON.NE.O)  GO  TO  10 
AM(2,1)=AM(1,2) 

AM(3 , 1)=AM( 1 ,3) 

AM ( 3 , 2 ) = AM (2,3) 

AM(4,1)=-AM(1,4) 

AM(4,2)=-AM(2,4) 

10  CONTINUE 

AS( 1 , 1)=S0RT( (AM(1 , 1)+AM(1 ,2)+AM(2 , 1)+AM(2, 2) ) /2. 0) 

AS  ( 1 ,  ?■ )  =SQRT  ( (  AM  ( 1 , 1 )  -  AM  (  i ,  2 )  +AM  (  2 , 1 )  -  AM  (  2 , 2 ) )  /  2 . 0 ) 

AS (2, 1;*SQRT((AM(1,1)+AM(1,2)-AM(2, l)-AM(2,2))/2.0) 

AS(2,2)=SQRT( (AM(1 , 1) -AM(1 ,2)-AM(2, 1)+AM(2 ,2) )/2. 0) 

C 

C  CALCULATE  THE  RELATIVE  PHASES  OF  (S) 

C 

PHI (1 ,2)=0 . 0 
X=AM( 1 , 3 )+AM(2 ,3 ) 

Y*AM(1,4)4-AM(2,4) 

PHI ( 1 , 1 )=PKI ( 1 , 2 )+ATAN2 ( Y , X) 

X=AM(3 , 1) -AM (3 ,2) 

Y=AM(4,1)-AM(4,2) 

PHI ( 2 , 2)=PHI ( 1 , 2 )+ATAN2 (Y , X) 

X=AM (1,3) -AM (2,3) 

Y»AMU,4)-AM(2,4) 

PHI(2,1)»PHI(2,2)+ATAN2(Y,X) 

C  PRINT, 'PKI(2,1)  *  ' , PHI (2,1) 

IF(MON.EQ.O)  PHI (2, 1)=PHI (1,2) 


ono  oooo  o  o  o  n  n 
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RETURN 

END 

C 

SUBROUTINE  AMMT0S(AM,AS,PHI) 

THIS  ROUTINE  CALCULATES  THE  RELATIVE  PHASE  SCATTERING  MATRIX  FROM 
THE  MODIFIED  MUELLER  MATRIX  IN  BISTATIC  AND  MONOSTATIC  CASES. 

PHI (1,2)  IS  ARBITRARY  E.G.  (  =  0.0  ) 

COMPLEX  Cl 

DIMENSION  AM(4,4),AS(2,2) ,PHI(2,2) 

COMMON  PI, PI 1, Cl, MON 
WRITE(6,*)PI,CI 

CALCULATES  THE  MAGNITUDES  OF  (S) 

IF (MON . NE . 0)  GO  TO  10 
AM ( 2 , 1 ) =AM (1,2) 

AM(3,2)=2.0*AM(2,3) 

AM(4,2)=-2.0*AM(2,4) 

10  CONTINUE 
DO  1  1=1,2 
DO  1  J=1 , 2 

AS(I , J)=SQRT(AM(I , J) ) 

1  CONTINUE 


CALCULATES  THE  RELATIVE  PHASES 


PHI(1,2)=0.0 
PHI ( 1 , 1)=PHI ( X , 

PHI ( 2 , 2 ) ^PHI ( 1 , 

IF (MON . EQ. 0)  GO 
Y=AM(3,4)+AM(4, 

X=AM(3 , 3) -AM(4 , 

PHI (2, 1)=PHI(1, 

GO  TO  21 

20  PHI (2, 1 ) — PHI ( 1 t 

21  CONTINUE 
RETURN 
END 

$ ENTRY 

0.500E  00-0. 494E  00  0.679E- 
0.491E  00-0 . 63 IE* 
0.500E  00-0. 403E  00  0.2.5  IE 
0.368E  00-0. 169E 
0.500E  00  0 .454E-01-0 . 606E- 
0.488E  00-C.815E- 
0 . 162E *02  0.471E-02  0.239E- 
0.983E  00  0.655E- 
0.306E-01  0 . 662E-01  0.413E- 
0.837E  00  O.UOE 
0.540E  00  0 , 573E-02’<0 . 344E* 
0.449E  00-0.262E- 


2)+ATAN2 (AM( 1 ,4) , AM( 1 ,3) ) 

2) +ATAN2(AM(4 ,2), AM (3, 2)) 
TO  20 

3) 

4) 

2)+ATAN2(Y,X) 


2) 


01-0. 
01  0, 
00-0, 
00  0, 
01  0, 
02  0. 
02  0, 
01-0, 
01  0, 
00-0, 
01  0 
01-0 


179E-01 
206E-01 
884E-01 
124E  00 
350E-03 
877E-01 
138E-02 
193E-0.1 
178E-01 
108E  00 
440E-01 
437E-Q1 
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1-3  THIRD  COMPUTER  PROGRAM  : 

This  program  reconstructs  the  scattering  matrix  [S]  from  the  knowledge  of  the 

spherical  coordinates  on  the  Poincare  sphere  of  two  C0P0L  nulls  or  one  C0P0L 

and  one  XPOL  null  using  the  equations  of  Section  3.7. 

//BAKRY  JOB 

/*J0BPARM  R=348 

//GO  EXEC  WATFIV , PARM=N0EXT 

//SYSIN  DD  * 

$J0B 

C 

C  THIS  PROGRAM  RECONSTRUCTS  THE  SCATTERING  MATRIX  (S)  WITH  RELATIVE 

C  PHASE  GIVEN  ITS  C0P0L  NULLS  (P, THETA, PHI ) ,  WHERE  P, THETA  AND  PHI 

C  ARE  THE  RADIUS,  COLATITUDE  AND  LONGITUDE  OF  THE  C0P0L  NULLS  ON  THE 
C  POINCARE'  SPHERE  OR  ONE  C0P0L  AND  ONE  XPOL  NULL. 

C  DATA:  PUT  NCC  (TWO  COPOL  NULLS)  FIRST. 

C  PUT  NCX  (ONE  COPOL  AND  ONE  XPOL  NULL)  SECOND. 

C 

DIMENSION  THETA(2) ,PHI (2) 

COMPLEX  Cl ,Q,RHO(2) ,R1 ,R2 , AK, S (2 , 2) , SI ,Q1 ,R3 ,C0NJG 
CI=(0. 0,1.0) 

PI=4.0*ATAH(1.0) 

PI180-PI/ 180.0 
PI 1=180 . 0/PI 
C  PRINT',  PI, Cl 
AMIN*l.E-03 
AMAX=1 .E+10 
NCC=12 
NCX=5 

NDATA=NCC+NCX 

DO  10  II-l.NDATA 

IF ( II . LE .NCC)  WRIT£(6 , 100) 

IF(II.GT.NCC)  WRITE(6, 110) 

READ  (5 , SO ) P , (fHETA(I) , PHI ( I ) ,1=1,2) 

WRITE ( 6 , 60 ) P , (THETA ( I ) , PHI (I) ,1=1,2) 

WRITE (6 , 70) 

DO  1  1=1,2 
THR=THETA(I)*PI180 
PHR=PHI(I)*PI180 
X=1 . O+COS (THR) 

Y=1.Q-C0S(THR) 

IF(Y.LE. 1 .E-07)  GO  TO  20 
Z=SQRT(X/Y) 

GO  TO  21 

20  Z=AMAX 

21  CONTINUE 

Q«Z*(C0S(PHR) -CI*SIN(PHR) ) 

Ql=1.0+Q 

C  PRINT, *  Q1  =  ' ,Q1 

IF (CABS (Ql) . LE. 1 .E-03)  GO  TO  30 
C  IF (CABS (qi).LE. AMIN)  Q1-1.0/AMAX 


ouuu 
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RHO(I)»-CI*(1.0-Q)/Q1 
GO  TO  31 

30  RHO(I)=-CI*AMAX 

31  CONTINUE 

WRITE(6,71)THETA(I) ,PHI(I) ,Q,RH0(I) 

1  CONTINUE 

IF(II .GT.NCC)  GO  TO  40 

C  PRINT, '**RHO(l) ,RH0(2)  =  ' ,RH0(1) ,RH0(2) 
Rl=RHO(l)+RHO(2) 

R2*RHO(l)*RHO(2) 

R1MAG=CABS(R1) 

X1=AIMAG(R1) 

X2=REAL(R1)  ^  _ 

IF  CABS (XI ) .LE. AMIN. AND. ABS(X2) .LE. AMIN)  GO  TO  45 
R1PHAS=ATAN2(X1,X2) 

GO  TO  46 

45  CONTINUE 
R1PHAS=0 . 0 

46  CONTINUE 

C  PRINT, '**TEST**' ,R1,R2,X1,X2,R1PHAS 
R2MAG=CABS(R2) 

X=RlMAG**2+2 . 0*R2MAG**2+2 . 0 

X=2 . 0*X 

AK=SQRT(P/X) 

S l=COS (R 1PHAS ) -C I*S IN (R 1PHAS ) 

S(1,1)=-2.0*R2*S1 

S(1,2)=R1MAG 

S(2,1)=S(1,2) 

S(2,2)=-2.0*S1 
GO  TO  41 
40  CONTINUE 


R}*Q(1) .  COPOL  NULL 

RH0(2) .  XPOL  NULL 


R l=RHO ( l ) **  2*C ON JG ( RHO ( 2 ) ) +RHO ( 2 ) 

R1MAG=CABS(R1) 

R1PHAS=ATAN2 (AIMAG (R1 ) ,REAL(R1 ) ) 

R2=RHO(1)-RHO(1)*CABS(RHO(2))**2-2.0*RHQ(2) 

R3=2 . 0*RH0 ( 1 )*C0N JG (RHO ( 2 ) ) -CABS (RHO ( 2 ) ) **2+1 . 0 
R2MAG=CABS(R2) 

R3MAG=CABS(R3)  _ 

D=2 . 0*R1MAG**2+(CABS (RHO (1) )*R2MAG)**2+R3MAG**2 
AK=SQRT(P/D) 

S1-*C0S(R1PHAS)  -CI*SIN(R1PHAS) 

S(l,  1)*RHQ(1)*S1*R2 
S(1,2)=R1MAG 
S(2, 1)*S(1 ,2) 

S(2,2)=-R3*S1 
41  CONTINUE 

WRITE (6, 80) AK 
DC  2  1*1,2 

WRITE (6 ,81) (S(I, J) , J-l ,2) 

2  CONTINUE 
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WRITE(6,90) 

DO  3  1=1,2 

WRITE (6, 81) (S(I, J)*AK,J=1,2) 

3  CONTINUE 
10  CONTINUE 
STOP 

50  F0RMAT(5F10.3) 

60  FORMAT  ('0' ,  'THE  INPUT  DATA  IS  :'//'o','P  =  '  ,F10.3/'0', 
l'THETA(l)  =  ' , F10 . 3 , 2X , ' DEGREES ' , 10X , ' PHI ( 1 )  =  ',F10.3,2X, 

2' DEGREES '/'O' . 'THETA(2)  '=  ' ,F10 . 3 ,2X , ' DEGREES ' , 10X, ’ PHI (2)  =  ' , 
3F10. 3, 2X, 'DEGREES'//) 

70  FORMATC'O' ,2X, ’THETACDEG.)' ,5X, 'PHI (DEG.) ' ,11X, 'Q' , 19X, 'RHO ' /) 

71  FORMATC'O' ,2(2X,F10. 3) ,2(2X,2E10 . 3) ) 

80  FORMATCl* ,  'THE  RECONSTRUCTED  SCATTERING  MATRIX  WITHOUT  MULTIPLYIN 
1G  BY  THE  CONSTANT  K  =  ' , 2F10 . 3 , 2X, ' IS  :  '//) 

81  FORMATC'O' ,2(2X,2F10. 3)) 

90  FORMATC'O' , 'THE  SCATTERING  MATRIX  IS  :  '//) 

100  FORMAT('l' , 'THE  COPOL  NULLS  ARE  KNOWN'//) 

110  FORMAT (' l' , 'ONE  COPOL  NULL  AND  ONE  XPOL  NULL  ARE  KNOWN'//) 

END 


SENTRY 

1.0 

109.414 

-4.055 

66.486 

-11.708 

1.0 

118.001 

-22.209 

41.311 

-45.001 

1.0 

172.57 

70.438 

7.736 

179.957 

1.0 

139.467 

-31.889 

40.877 

58.378 

1.0 

177.921 

-28.335 

0.934 

-31.12 

2.0 

0.0 

180.0 

2.0 

90.0 

-90.0 

90.0 

90.0 

1.0 

180.0 

180.0 

1.0 

0.0 

0.0 

1.0 

90.0 

180.0 

90.0 

180.0 

1.0 

90.0 

0.0 

90.0 

0.0 

1.0 

90.0 

270.0 

90.0 

270.0 

1.0 

109.414 

-4.055 

92.054 

172.172 

1.0 

118.001 

-22.209 

100.543 

148.062 

1.0 

172.57 

70.438 

89.735 

-53.177 

1.0 

139.467 

-31.889 

89.756 

-166.555 

1.0 

177.921 

-28.335 

90.568 

150.796 

DISTRIBUTION  LIST  (A) 


Commander 

Air  Force  Avionics  Laboratory 
ATTN:  Mr.  Allen  Blume  (AFWAL/AADM) 

Mr.  Robert  L.  Davis  (AFAL/WRP) 

Mr.  Harold  Weber  (AFWAL/AADM) 
Wright-Patterson  Air  Force  Base,  OH  45433 

ETE 

ATTN:  Dr.  A.  Schell  (RADC/ETE) 

Hanscome  Air  Force  Base,  MA  01731 

General  Electric/RSD 

ATTN:  Mr.  A.  B.  Graf inger/Room  6258H 

3193  Chesnut  Street 

Philadephia,  PA  19101 

R.C.  Hansen,  Inc. 

Box  215 

Tarzana,  GA  91356 
ITT  Gilfillan 

ATTN:  Dr.  David  E.  Hammers 
7821  Orion  Avenue 
P.0.  Box  7713 
Van  Nuys,  CA  91409 

Commander 

Naval  Avionics  Center 
ATTN:  Mr.  Paul  Brink 
21st  &  Arlington  Avenue 
Indianapolis,  IN  46218 

Superintendent 

Naval  Postgraduate  School 

ATTN:  Dr.  Lonnie  Wilson  (Code  62W1) 

Monterey,  CA  93940 

Commanding  Officer 

Naval  Research  Laboratory 

ATTN:  Mr.  John  Daley  (Code  7946) 

Mr.  Fred  Staudaher  (Code  5368) 

Mr.  Denzil  Stilwell  (Code  7910) 

Commander 

Naval  Sea  Systems  Command 
ATTN:  Mr.  C.  Jedrey  (SEA-62R13) 

Washington,  DC  20362 


DISTRIBUTION  LIST  (A)  PAGE  2 


Commander 

Naval  Surface  Weapon  Center 
ATTN:  Dr.  B.  Hollmann  (Code  F12) 

Mr.  John  Teti  (Code  DF34) 
Dahlgren,  VA  22448 

Commander 

Naval  Weapons  Center 
Electronic  Warfare  Department 
Microwave  Development  Division 
ATTN:  Mr.  F.  F.  St.  George  (Code  354) 
China  Lake,  CA  93555 

Office  of  Naval  Research 
ATTN:  Dr.  H.  Mullaney  (Code  427) 

800  North  Quincy  Street 
Arlington,  VA  22217 

The  Ohio  State  University 
ElectroScience  Laboratory 
ATTN:  John  Young 
1320  Kinnear  Road 
Columbus,  OH  43212 

Commanding  Officer 

Rome  Air  Development  Center 

ATTN:  Mr.  Thomas  Shields  (OCDR) 

Griffis  Air  Force  Base,  NY  13440 

Technology  Service  Corporation 
ATTN:  Dr.  Fred  Nathanson 
8555  16th  Street,  Suite  300 
Silver  Spring,  MD  20910 

Teledyne  Micronetics 
ATTN:  Dr.  Steven  Weisbrod 
7155  Mission  Gorge  Road 
San  Diego,  CA  92120 

Commander 

US  Army  Electronics  Command 
CS  and  TA  Laboratory 
ATTN:  Mr.  Boaz  Gelerneter  (DRSEL-CT-R) 
Mr.  Willy  Johnson 
Ft.  Monmouth,  NJ  07703 

Commander 

US  Army  Missile  Command 

ATTN:  Dr.  James  Wright  (DRSMI-RER) 

Redstone  Arsenal,  AL  35809 


DISTRIBUTION  list  (b) 


Defense  Documentation  Center 
Cameron  Station 
Alexandria,  VA  22314 

Director 

Naval  Research  Laboratory 
4555  Overlook  Avenue,  SW 
ATTN:  Tech.  Information  Officar 
Washington,  DC  20375 

Director 

Naval  Research  Laboratory 
4555  Overlook  Avenue,  SW 
ATTN:  Library,  Code  2620 
Washington,  DC  20375 

Office  of  Naval  Research 
Mathematics  Program,  Code  432 
Arlington,  VA  22217 

Office  of  Naval  Research 
New  York  Area  Office 
715  Broadway-5th  Floor 
New  York,  DY  10003 

Command  Officer 

Office  of  Naval  Research/Branch  Office 
536  South  Clark  Street 
Chicago,  I'L  60605 

Command  Officer 

Office  of  Naval  Research.  Eastern/ 
Central  kegional  Office 
Bldg  114,  Section  D 
666  Summer  Street 
Boston,  MA  02210 

Command  Officer 

Offjce  of  Naval  Research  Western 
Regional  Office 
1030  East  Creen  Street 
Pasadena,,  CA  91106 


DISTRIBUTION  LIST  (B)  PAGE 


Professor  W.  Boerner 
University  of  Illinois 
N00014-80-C-0708 
NR  041-565 

Dr.  J.  DeSanto 

Naval  Research  Laboratory 

Code  4340 

Washington,  DC  20375 

Dr.  J.  Heacock 

Office  of  Naval  Research 

Code  463 

Arlington,  VA  22217 
Dr.  G.  Heiche 

Naval  Air  Systems  Command  Headquarters 
Jefferson  Plaza  #1 
1411  Jefferson  Davis  Highway 
Arlington,  VA  22516 

Dr.  A.  Jordan 

Naval  Research  Laboratory 

Code  7943 

Washington,  DC  20375 

Dr.  M.  McKisic 

Office  of  Naval  Research 

Code  486 

Arlington,  VA  22217 

Dr.  H.  Mullaney 
Office  of  Naval  Research 
Code  427 

Arlington,  VA  22217 
Mr.  J.  Willis 

Naval  Air  Systems  Command  Headquarters 
Jefferson  Plaza  1 
1411  Jefferson  Davis  Highway 
Arlington,  VA  22516 

Dr.  Walter  K.  Kahn,  Professor 
Department  Chairman 

Electrical  Engineering  and  Computer  Science 
School  of  Engineering  and  Applied  Science 
The  George  Washington  University 
Washington,  DC  20052 


DISTRIBUTION  LIST  (C)  DOMESTIC 


Dr.  Robert  G.  Appling 
Senior  Research  Engineer,  RAIL 
Engineering  Experiment  Station 
Georgia  Institute  of  Technology 

Dr.  David  Atlas  (Code  910) 

Laboratory  for  Atmospheric  Sciences 
NASA,  Goddard  Space  Flight  Center 
Greenbelt,  Maryland  20771 

Dr.  Ivan  Cindrich 
Radar  and  Optic  Division 
ERIM,  3300  Plymouth  Rd. 

P-0.  Box  8618 

Ann  Arbor,  MI  48107 

Dr.  Dick  Doviak 

Chief,  Advance  Techniques 

National  Severe  Storms  Laboratory 

1313  Halley  Circle 

Norman,  OK  '73069 

Dr.  Adrian  K.  Fung 

Department  of  Electrical  Engineering 
University  of  Kansas 
Lawrence,  KS  66045 

Dr.  Edmond  Gillespie, 

Chairman 

Electrical  and  Computer  Engineering  Dept. 
California  State  University  Northridge 
18111  Nordhoff  Street 
Northridge,  CA  91330 

Dr.  Marlin  Gillette 
Mailstop  B-5C 
Bell  Aerospace  Textron 
P.0.  Box  1 
Buffalo,  NY  14240 

Dr.  J.W.  Goodman 

Department  of  Electrical  Engineering 
Stanford  University 
Stanford,  CA  94305 

Dr.  Doran  Hess,  Chief  Physicist 
Scientific  Atlanta 
3845  Pleasantdale  Road 
Atlanta,  GA  30340 


DISTRIBUTION  LIST  (C>  DOMESTIC 


Dr.  J.  Richard  Huynen 

10531  Blandor  Way 

Los  Altos  Hills,  CA  94022 

Professor  Dr.  Akira  Ishimaru 
Department  of  Electrical  Engineering 
University  of  Washington 
P.0.  Box  88062 
Seattle,  WA  98188 

Dr.  Edward  M.  Kennaugh 
ElectroScience  Laboratory 
Department  of  Electrical  Engineering 
The  Ohio  State  University 
1320  Kinnear  Road 
Columbus,  OH  43212 

Dr.  Wolfgang  H.  Rummer 
1310  Sunset  Avenue 
Santa  Monica,  CA  90405 

Dr.  Louis  N.  Medgyesi-Mitschang 
Senior  Scientist 

McDonnell  Douglas  Research  Laboratory 
Box  516 

St.  Louis,  M0  63166 

Dr.  James  Metcalf 
Radar  Applications  Group 
Georgia  Institute  of  Technology 
Electrical  Engineering  Department 
Atlanta,  GA  30332 

Dr.  E.  K.  Miller 

Electronics  Engineering  Department 
University  of  California 
Lawrence  Livermore  Laboratory 
P.0.  Box  808 
Livermore,  CA  94550 

Dr.  Charles  Miller,  Director 
Electromagnetics  Division 
National  Bureau  of  Standards 
Boulder,  CO  80303 

Dr.  Allen  Taflove 
Research  Engineer 

EM  Interaction  with  Complex  Systems 
IIT  Research  Institute 
10  West  35th  St. 

Chicago,  IL  60616 


DISTRIBUTION  LIST  (C)  DOMESTIC 


Dr.  Gus  Tricoles 
General  Dynamics 
Electronics  Division 
P.O.  Box  81127 
San  Diego,  CA  92138 

Dr.  Charles  Warner 

Department  of  Environmental  Sciences 

Clark  Hall 

University  of  Virginia 
Charlottesville,  VA  22903 


